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Chapter  I 


INTRODUCTION 

Inuadatlcn  oD  river  flood  plains  by  occasional  flood 
discharges  has  been  a  threat  to  the  life  and  propei’ty  of  Kan 
since  the  davjD,  of  civilization^  Nan  has  long  sought  maars 
of  predicting  floods  and  averting  thair  clanger*  liodern  man’s 
need  for  quantitative  information  concerning  flood  flous 
becomes  even  more  pressing  as  he  alters  natural  waterways 
vrith  structUTes,  realignments,  and  diversions* 

A  ms-Jor  problem  continually  faced  by  highway  designers,  urban 
planners,  and  watershed  management  engineers  is  that  of  ceteroilning  the 
frequency  of  peak  flood  discharges  from  stream  and  river  basins*  For 
large  river  basins  (in  excess  of  100  square  miles)  flood  frequency  in- 
foriaation  is  generally  sivailable  by  virtue  of  long  term  records  of  pre¬ 
cipitation,  stream  flov?,  etc*,  which  have  been  collected  by  various 
public  agencies  (includixjg  the  U*  S,  Environmental  Science  Services 
Administration,  Geological  Survey,  Bureau  of  Reclamation,  Agriculture 
Research  Service,  and  others)*  For  watersheds  smaller  than  100  square 
iniles  ty.i  area  there  is  generally  a  shortage  of  hydrologic  data*  Only 
within  the  jjast  few  years  has  significant  research,  been  dii’ectcd  toward 
a  study  of  flood  frequencies  on  basins  of  this  size. 

The  ‘’D2:'aij:age  Correlation  Research  Froject",  which  was  initiated  by 
the  Department  of  Civil  Engineering  and  Engineering  Mechanics,  Montana 
State  lliiversit^r ,  Bozeman,  Montana,  in  1963;  addresses  itself  to  the 
problera  of  predicting  the  frequency  with  which  a  peak  discharge  of  given 
magnitude  may  be  expected  on  sirall  watersheds  In  Montana*  This  invesb?.- 
gat  ion  is  sponsc>red  by  the  Montana  State  Highway  Department  and  the  U*  S* 
Bureau  of  Public  Roads.  The  work  is  being  done  under  the  direction  of 


7T. 


A  ...  J 

4.  I  -  ■ 


t  * 

*1  ‘‘• 

\ 


•  -v' 


>  ♦*“  • 

'  I  IM 


«. 


>  e4*  if-.  ■' \t 


i 


*f  I'l  "'O  *.^v'.  •■  •:  '  "^v'*  ■' 

%  •.  -  ;v.  . . 

j*..'  ’  i’’  '.  :t<*’ ‘.'  ' 

lk‘  \/'r  1  ■ 

-ii '  (- '  ^  •  _!,  •  I  ■•  li  * 

^  X  ■ .  ■)  '  "iA  n  ?  ■  *  •■•.♦ 


«*  ' 


>•<  '%  :•!/. 


.y  . 

?'  .  vrfV  '■  ' 

^  ■.  ^  ■  t  ;  'lf> 

'-  ■'■■  ;  54i!i 

.  .■;  :  :4 

.-I 


■^<  •.  /'’ 


-  1  u 

i' 


■'  .f  yr-f 


•1 


%  i 

\  y  'ijc  .  /u  N’  ‘-/'v  ^•‘Wj^i  f  A 


cj..  .  .4  -•  ».*■  ri.'  '.;.  .  -■>/'•  I '.p  i*  ^rv^i  fim^k'  '♦i.Cq 

*  f, 

i'*  4  .'’k:  '  >'i’.  ff  '^O- '•  *f.  H. 


■"« 


■'  ,«  d  iWS  ■'*’'■  ’^'^  .■-•'•A  Ui)  ^11*4,4 

e  M  f  v-  t  '  ,  •  ■•  I  •  ,'rf,  c  r  '*^  :  ’  V  ftl  'it)l 

’^4  i  n  if.  i'S  ,  ,,TLi;g 


2  •.  'i V-"  'r  . ■  *.  .'I.  .4,^13  .'\#. ->>) 


<  .*  '  I 


M 


')J  p'.^  •  .'/ -C  U‘.H'-(.!U  ^',  .  >/ '.i^^  “|.,t>  M 


*17  *  ■  .  , 

'  '  '.  “J  •xfil.t''  “  .  »tfv  if\ii  V'  ■-r.  'T  .  '  ‘  rij  '.i’ij'i  '/Tc?!’  ({s  '  '  '.■>»  '^H.' 

»  ■  ■  M  ,T  - «  ’  . 

i  V 

’CSi^t)  '^'ii  ‘K*  O'Ib  ’.  cl-^“ •/•i'  vHXi'a 

'■  U  *'■  y 

>i.0'«X4  >  tii  '  X  v' -;  rv^‘.  2  tid.  ’j)  ~  i  k*i 

y  *  . 


*  • 


®-r 

■t 


'44*  ../•''I..  V  ^-'-^-I^  S,  vtr'l 


4  • »-•</  '"L'-  -  i.t  ^‘j-^-'i  '.}.  •'  iX'ftV'i  ' '  \  • 

k  . 

'i  1  :  -:-tBX.<i  ;i^i>*.;f  V  y.j  nf  a^'^w  • 

*tovi/l  4?i!fr  .  -•nA!  •’!  nt  -  **.  i 

.If  ^4.ii  ti^A  Jffsw^ ‘,y\  ,  4  • 

■  .  *'*:  ■  Vj5^'. i''-^  • .'  ^>’J~  ,-i  ■ 


4 

fn 


L^-t  I 


.1 


^  ti* 

W-Jr' t  ^ 

,  litt--'.'  > 

'"4  .i'-'x  1 

•  Jl  V'i 

*  ■>  ‘ 

*  X 

JWH  ^ 

iX^  •  tp-Xdo'f<i 

..  i  i*  ii, 

•'A  /■■'*/ 

3^-(r  •  '3 

bJi 

•ifi  V  ?f  tNtg:-: 

Jk 

•vr.^kl  .  i; 

■»5p^  ’V'  .1  - 
*  • 

1'  5t' ^'3 

>3', ';*r''  ?I 

yijXSi.i  3^; 

J 


li'.  f . 

i'-l4,..  , 


-  6  - 


Theodore  T,  Williams,  Associate  Professor  of  Civil  Engineering  and  En¬ 
gineering  Machanics  at  Montana  State  University. 

The  Drainage  Correlation  Research  Project  is  a  two  phase  study: 
a)  to  determine  flood  frequency  relationships  on  watersheds  smaller  than 
100  square  miles  from  long-term  climatological  data  already  in  existence; 
and  b)  to  make  comprehensive  hydrologic  studies  of  four  small  watersheds 
in  eastern  Montana,  The  watersheds  selected  for  study  are  Bacon  Creek 
in  Wheatland  County,  Duck  Creek  in  Prairie  and  McCone  Counties,  Hump 
Creek  in  Sweet  Grass  County,  and  Lone  ^hln  Coulee  in  Pondera  County,  Con¬ 
tinuous  records  of  climatic  factors  and  streamflou  are  being  obtained  from 
these  four  basins.  A  variety  of  peak-flow  prediction  techniques,  which 
utilize  data  such  as  that  being  collected,  are  being  examined. 

Among  the  techniques  being  considered  under  the  second  phase  of  the 
Drainage  Correlation  Research  Project  is  one  which  was  developed  by  the 
Soil  Conservation  Service,  U.  S.  Department  of  Agriculture  (SCS).  The 
SCS  method  posse'sses  mny  desirable  features,  and  seems  to  have  consid¬ 
erable  potential  as  a  tool  for  predicting  flood  frequencies.  The  study 
reported  in  this  report  analyzes  the  method  developed  by  the  SCS;  describes 
the  use  of  a  related  computer  program;  tests  the  program's  ability  to  re¬ 
produce  actual  hydrographs  from  a  selected  watershed;  and  evaluates  the 
possible  utility  of  this  method  on  Montana  watersheds. 

The  investigation  reported  herein  consisted  of  the  analysis,  using 
the  SCS  method,  of  a  single  rain-caused  runoff  event  which  occurred  on 
Duck  Creek,  the  largest  of  the  four  watersheds  being  studied  by  the 
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I)ra?lnage  Correlo-tion  Eesearch  Project «=  Pack  Creok  is  an  ephsLiaral  stream 
typical  of  most  east  cm  Montana  streams  which  drain  sicall  watersheds*  It 
is  dr^’-  most  of  the  time^  but  occasionally  has  surface  flow  after  <a  raiii” 
caused  runoff  eveiiu  or  canning  spring  snowmelt.  In  the  four  years  for 
which  data  have  been  obtained  at  Duck  Creek;?  there  have  been  only  a  few 
runoff  events  of  consequence^  and  only  an  evenij  which  occurred  on  June  ?.6, 
1965  lends  itself  tc  an  analysis  by  the  SC3  methed* 

A  nujaber  of  events  v?Quld  need  to  be  analyzed  before  conclusive  state- 
inants  as  to  the  validity  of  the  SCS  method  could  be  made*  IIaverthelesG_? 
it  is  believed  that  the  results  of  the  analysis  of  this  one  storm,  will  be 
a  valuable  contribution  in  estimating  the  ai^plicabllity  of  the  iriethod  to 
Mont  a  na  w’at  0  rshe  d  s  * 


Chapter  II 


LITERATURE  REVIEW 

Scientific  hydrology  is  a  relatively  new  area  of  study.  Appli¬ 
cation  of  the  scientific  approach  to  hydrologic  problems  has  had  its 
greatest  growth  in  this  century,  and  especially  in  the  last  thirty 
years . 

Hydrologic  measurements  were  not  considered  to  be  of  much  import¬ 
ance  during  the  early  history  of  the  United  States.  Population  centers 
developed  along  the  rivers  and  lakes  and  so  no  shortage  of  water  was 
experienced  until  the  westward  expansion  into  the  arid  and  semi-arid 
regions.  The  beginning  of  systematic  collection  of  hydrologic  data  in 
this  country  can  probably  be  taken  to  be  in  l888,  when  the  U.  S.  Geo¬ 
logical  Survey  under  the  direction  of  Frederick  H.  Newell  set  up  its 
first  river- measurement  station  on  the  Rio  Grande  R;Lver  at  Embudo,  New 
Mexico. 

Modern  Hydrology 

Since  about  1930,  the  increase  in  the  amount  of  printed  information 
made  available  in  the  field  of  hydrology  has  increased  very  rapidly. 

This  indicates  something  of  the  increased  need  there  has  been  in  recent 
years  to  obtain  more  data  and  develop  better  prediction  methods. 

The  study  of  hydrology  can  be  subdivided  into  the  two  general  areas 
of  stochastic  and  deterministic  hydrology.  Deterministic  hydrology  can 
be  further  separated  into  physical  or  analytical,  dynamic  and  parametric. 

Stochastic  hydrology  attempts  to  predict  future  events  or  reconstruct 
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past  events  based  on  the  statistical  properties  of  the  knovn  record. 

White  {19^7 )  states  that  if  the  value  that  a  variable  has  implies  an 
element  of  chance^  then  it  is  a  stochastic  varirablsc  Thus  if  the  2‘ecord 
aval.lablo  for  some  variable  (the  value  of  which  is  random  in  nature  or 
stochastic)  is  restricted  in  tinie^  then  the  predicted  futuxe  projection 
of  this  i’ecord  is  done  by  statistical  means » 

Physical  (amlytical )  hydrology  addresses  Itsa.li'  to  particu,lar 
specialised  •probloms  of.  the  hydrologic  process.  It  does  not  try  to  supply 
any  ansv^ers  to  peak  discharge,  annuaj.  yield,  or  frequency  slnidicGe  In 
a  sense  this  branch  of  hydrology  is  pure  research  in  that  it  seeks  to 
establish  relationships  betu’een  variables  operating  in  the  hydrologic 
cycle  but  does  not  attempt  to  soj.ve  any  related  practical  prob.lemx. 

Pynani.ic  hydrology  is  the  term  for  hiydrologic  studies  involving  the 
dynamic  v?ave  theories  of  fluid  flow.  These  wave  theories  have  been 
applied  to  overland  flow  as  well  as  open  channel  flow  to  derive  syn“ 
thctic  runoff  byclrographs. 

Parainetric  hydrology  attempts  to  discover  the  relat ionsbips  among 
physical  jjarameters  that  are  involved  in  the  particular  hydrolcgic  events 
that  are  of  interest  and  than  to  \i3e  them  to  simulate  non^ recorded  events. 
Amorocho  and  Haa-t  (1964)  list  methods  of  correlation  analysis,  partial 
system  synthesis  with  linear  or  nonlinear  analysis,  and  general  system 
synthesis  as  metheds  used  in  paraBietric  hydrology.. 

In  the  prediction  of  peo.k  runoff  fi'om  a  small  watershed,  by  a  method 
such  as  that  developed  by  the  SCS,  both  stochastic  and  deterministic 
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hydrology  come  into  play.  The  streamflow  itself,  being  the  result 
of  precipitation,  is  a  stochastic  process.  Watershed  characteristics, 
vhich  modify  or  affect  the  streamflow,  are,  by  themselves,  physical 
parameters.  Dynamic  hydrology  must  be  considered  in  the  flood  routing 
process  while  parametric  hydrology  is  utilized  to  find  correlations 
among  the  various  watershed  parameters, 

A  discussion  of  the  technique  for  predicting  flood  frequencies 
which  is  currently  in  use  by  the  Soil  Conservation  Service,  is  reserved 
for  Chapter  III,  because  it  is  the  basis  for  this  report. 


Chapter  13,1 


SOIL  COLSLRVATIOI  SEBVICS  RUIIO?F  PREDICTION  METHOD 
The  Soil  Cons er'/at ion  Service  method  of  predicting  runoff  from  im** 
gaged  v?atersheds  is  characterized  by  the  development  of  a  synthetic  unit 
hydrography  Under  this  raethod  the  storm  characteristics  are  transformed 
into  a  synthesized  f3,ood  runoff  hydrograph  by  the  basin  characteristics^ 
the  shape  of  the  synthetic  unit  byclrograph  and  the  baseflov/  hydrography 
The  general  procodure  which  is  followed  in  developing  a  flood  runoff 
h;ydrograph  is  shown  by  the  block  diagram  in  Figure  Ic 

As  showui  in  Figure  the  storm  runoff  voluna  must  first  be  derived 

from  the  storm  characteristics  and  the  basin  characteristics o  This  storm 
runoff  volume  is  then  shaped  into  a  storm  hydrograph  by  use  cf  a.  Byn- 
thetic  unit  hydrograph  whi.ch  in  turn  is  derived  from  the  basin  and  storiri 
characteristics*  A  bascfi.ow  hydrograph  then  is  added  to  the  storm  hy- 
drograph  as  indicated  in  Figure  1  to  produce  the  f3.cod  runoff  hydrograph* 
This  sequence  of  operations  involved  in  producing  the  flood  runoff  by- 
drograph  comprises  a  irmthematical  watershed  model*  Tiie  general  frame- 
W’ork  of  the  \fatershed  model  has  been  incorporated  into  a  program  for 
solution  by  a  digital  com'puter*  The  basin  characteristics  of  a  parti- 
cu3,ar  v/aterslisd  must  be  supplied  as  input  to  the  computer  to  transform 
the  general  model  into  a  specific  model  for  the  watershed  in  quest ion* 
FTOT)Osed.  Use  SCS  Ih'th^^ 

Tile  propOvSed  way  of  testing  the  ECS  method"’  is  to  use  it  on  an  actual 


'"This  method  is  outlined  In  the  iU'^’^i^rology  section  of  the  SCS  Rational 
Eng’lneering  llandbcok  (.1964)* 
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storm  where  the  input  (storm)  and  the  output  (runoff  hydrograph)  are 
known.  The  procedure  should  be  to  supply  data  characterizing  the  water¬ 
shed  and  storm  characteristics  as  input  to  the  computer  program  which  then 
uses  the  watershed  model  to  simulate  a  runoff  hydrograph.  The  volume  of 
the  simulated  hydrograph  is  compared  with  that  of  the  actual  event.  If 
there  are  major  discrepancies  in  the  volumes,  the  values  characterizing 
the  basin  parameters  may  be  altered,  and  the  computer  program  run  again. 
V/hen  the  volumes  are  adjusted  to  be  essentially  the  same,  the  model  is 
considered  to  be  a  valid  representation  of  the  watershed.  The  validity 
of  the  SCS  method  depends  upon  the  ability  of  the  adjusted  model  to 
accurately  reproduce  the  actual  runoff  hydrograph,  including  the  peak 
discharge. 

After  the  model  has  been  adjusted,  it  is  then  possible  to  route 
design  storms  through  the  model  and  thereby  simulate  design  hydrographs. 
Unit  Hydrograph  Theory 

The  unit  hydrograph  (UH)  for  a  given  watersVied  is  defined  as  the  dis¬ 
charge -to-time  relationship  that  yields  one  inch  of  runoff  from  a  storm 
of  given  duration  over  the  entire  watershed  area.  Unit  hydrographs  are 
generally  derived  from  as  many  actual  recorded  storms  as  possible  and 
then  generally  are  assumed  to  represent  the  unit  runoff  functions  for 
all  storms  of  similar  durations.  Sherman  (1932)  is  given  credit  for  for¬ 
mulation  of  the  unit  hydrograph  theory  while  Snyder  (I938)  first  intro¬ 
duced  a  method  for  constructing  a  synthetic  UH  which  may  be  used  for  the 
study  of  ungaged  watersheds.  Several  other  methods  for  defining  synthetic 
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UH  have  been  d,ovc3.cpsd  since  Snyder  iiiclu/Ilns  tliosc  by  ComonG 
Mlt-eliell  (19-’:'8),  and  Gray  (I961), 

System  Syntliesis  ayid  Anadys ts 

d/b,ere  no  actual  hydrographs  arc  available  for  a  vatcrshscl,  the  SCS 
method  of  i^rodicting  I’linoff  can  be  claL-eified  as  a  method  of  general 
system  synthesis ^  If  an  actual  unit  hydrograph  is  availabJ^o,  on  the 
other  hand;,  the  SCri  rathcd  v/ould  be  one  of  partial  system  sjynthosis  uith 
linear  system  analysis  * 

System  synthesis  accoj’ding  to  Amorocho  and  Hart  (1964)  Is  a  msthoa 
of  describing  the  operation  of  a  physical  system  uith  a  combination  of 
components  that  exist  in  the  system  and  uhose  functions  are  knonui  and  pre< 
dictable* 

System  ana3.ysis,  on  the  other  hand^  is  a  method  by  uhich  the  re¬ 
lationship  between  the  input  and  output  to  the  system  is  established 
Tiathemat ically  by  measui'iug  only  the  properties  of  the  input  (storm  data) 
and  output  (runoff)  v/ithout  regarding  the  nature  of  the  syGtern» 

An  exp-lanatioii  of  the  linear  and  nonlinear  properties  of  runoff 
prediction  raethods  will  follow  later  in  this  chapter. 

On  a  gaged  watershed  since  an  actual  UH  is  available^j  it  is  an 
anaJ.ytic  (as  opposed  to  synthetic)  function  since  the  internal  charac¬ 
teristics  of  the  system  are  not  known  or  synthesizedo  Even  though  the 
actual  hydrograph  is  available,  however,  synthetic  modifications  are 
necessary*  Interception  of  preciimtation  by  vegetation  and  baseflow 
characteristics,  for  instance,  must  be  synthesized*  Therefore,  even  on 
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a  gaged  \:atershGd^  the  GCS  method  is  partially  synihetic^  and  may  be 
classified  as  one  of  partial  system  synthesis  vrifh  linear  ana.lysis 


On  an  ungaged  V/atershed^  uhere  It  is  necessary  to  use  a  synthetic 
UH  the  entire  system  is  synthetic;  t'mxs  the  classification  of  general 
synthesis  method  applies  to  the  prediction  of  runoff  hydrographs  on  un- 
gaged  watersheds  £ 

Funct ioning  of  a  General  Synthesis  t!athod 

The  system  which  transforms  rainfall  to  runoff  can  be  thought  of  as 
made  uj)  of  three  separate  subsystems  when  the  lurit  hydrograph  approach 
is  being  considered*  •  1)  Ihe  first  subsystem;^  vhiich  creates  a  ’’rainfall 
excess"  function,  modifies  the  total  rainfall  input  to  account  for  in- 
f i.ltration,  interception,  and  depression  storage*  The  rainfall  excess 
thus  computed  is  the  amount-  of  rainfall  that  is  available  for  runoff, 

2)  The  second  subsystem  creates  a  "storm  runoff"  fimction  by  operating  on 
the  rainfall  excess  function.  This  subsystem  Teakes  use  of  the  basin’s 
topographic  characteristics,  3)  Finally,  the  third  subsystem,  which  cre¬ 
ates  the  "flood  runoff"  function  uses  baseflov?  infoimoation  to  alter  the 
storm  runoff  funct ion.  The  combined  effect  of  these  three  subsysteiss  is 
assum.cd  to  duplicate  the  natural  processes  which  occur  on  the  actual 
v'atershcd  being  studied.  The  paragraphs  that  follow  describe  in  detail 
these  three  subsystems  as  synthesized  by  the  SCS  method, 

FalnfalX  Excess  i  Tlie  ramfall  excess  function  (x^rc-duced  by  the  first 
subsystem)  is  given,  according  to  the  SCS  method,  by  equation  (3)» 
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(3) 


v/hare  Q  is  the  aiiount  of  rainfall  excess  in  inches  over  the  v.’atarshed_, 
P  is  the  storm  rainfall  in  inches,  and  S  is  defined  as  the  maximum 
potential  difference  between  P  and  Q  (hence  the  inaxiinum  Infiltration, 
capacity)  at  the  time  of  the  storin' s  beginning. 

Equation  (3)  Is  cased  on  a  hypothesis.  If  the  equD.tion  can  be 
shown  to  be  true,  then  the  hypothesis  can  be  assumed  to  be  true.  The 
original  hypothesis  can  be  stated  as  such: 


—  and  Q*»>1  as  P 
S  P 


(^0 


where  G  is  tbs  actual  retention  during  a  storm,  S  is  the  potential 
maximum  retention,  Q  is  the  direct  runoff  (or  the  actual  runoff),  and  P 
is  the  total  storm  raiirfall  (or  the  potential  maximum  runoff).  So  when 
flood  producing  storms  are  considered,  it  can  be  said  that: 


G 

S 


Q 

P 


(U 


Equation  (5)  cannot  be  proven  mathematically,  however,  it  folloi.’s  from 
equation  (3)  which  empirically  has  been  found  to  be  valid. 

Since  G  =  P  -=  Q,  coiuition  (5)  can  be  rewritten: 


P  -  Q 
S 


Q 

P 


(6) 


The  SC3  K^'drolcgy  randbooh  begins  the  devoJopment  cf  equation  (3)  with 
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tbc  Btatemont  of  equation  (6)6 

As  P  goer,  t-o  Infinity,  the  ratio  (P  Q)/S-^'l  and  oyp->l,  Keei')- 
ing  in  mind  that  P,  Q,  and  S  are  total  volumes  for  a  storm,  (P  “  Q) 
finally  fills  the  entire  soil  profile  to  its  maximum,  value  S  as  de¬ 
termined  for  the  condition  of  the  watershed  at  the  beginning  of  a  stoiTii 
Nov.’,  solving  directly  for  runoff  volu’ne  from  equation  (6); 


Q  ^ 


p  c 


(7) 


IJ 


A  further  adjustment  of  P  can  be  made  for  the  Initial  abs bract  ion 
Iq  which  reduces  P  by  the  amount  of  Interception,  depression  stcn.’age, 
and  infiltration  at  the  storm’s  beginning  before  runoff  occurs.  Several 
relationships  can  noi^  be  revised  to  a3.l0/.’  for  the  initial  abstraction 
(in  effect,  this  amounts  to  redefining  terms  v.’hieh  nov;  have  slightly 
diffci’ent  meanings  than  the  corresponding  terms  in  equations  (b)  through 

(7)  ). 


r  ::  Q  G  la 


(8) 


S  -i  G  -  P  -  Q 


(9) 


(P  -  la)  Q 


(10) 


Equation  (7)  can  also  be  restated 


P  »  la-f  S 


(11) 


Existing  data  from  several  watersheds  V7ere  studied  and  it  v.’as  found  chat 


N 
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can  be  taken  to  be  equal  to  .25 »  IBy  cubstltutin,"  this  va3j.ie  in 
equation  (11 equation  (3)  is  finally  derived. 

The  SCS  has  prepared  a  graijliical  solution  (runoff  as  a  function  of 
precipitation)  of  equation  (3)^  v/hich  is  shoun  in  Appendix  A.  As  is 
evident  in  Appendix  A,  precipitation  is  related  to  runoff  through  a 
fami.ly  of  curves,  each  curve  being  dram  for  a  particular  curve  number 
(CN).  'ilie  curve  numbers  are  a  function  of  G  and  are  given  by  the  re*= 
lat  tons  hip ; 


CN  r. 


1000 


Q  .L 
u  r 


10 


02) 


The  pararueter  CIl  is  used  instead  of  S  as  integer  valu.cs  of  S  do 
not  produce  curves  that  lend  themselves  to  easy  interpolation.  O'.he 
curve  numbers,  hov7ever,  make  the  P  vs.  Q  plot  easier  to  use. 

The  runoff  equixtion,  equation  (3),  used  by  the  SCS  is  a  relationship 
betueen  P  and  Q  involving  only  one  paranveter,  the  parameter  Vjaing  S.  S 
can  non  ba  empirically  related  to  as  iisany  characteristics  as  is  thought 
to  be  necessary.  The  developers  of  this  method  found  that  soil  and  cover 
conditions  had  the  most  effect  on  the  value  of  S  (and  therefore  on  the 
Cil)®  fhe  curve  number  CN  is  therefore  termed  the  soil  cover  corirjlex 
number . 


For  various  combinations  of  soil  type  and  vegetal  cover,  soil  cover 
complex  numbers  have  been  developed  cmj'iirically.  To  determine  such  a 
number,  srnll  uatersheds  Vi’ere  found  v’hich  had  onJ.y  one  typo  of  soil  a7id 
only  one  type  of  cover  condition,  A  number  of  points  were  plotted,  one 
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point  par  store:,  on  e  graph  of  total  precipitation  in  inchss  vc-rsus  total 
runoff  in  inches,  Tha  Oil  curve  that  best  fits  this  plot  is  then  taken  to 
be  the  representative  enrve  number  for  that  soil  cover  complex.  The.  re¬ 
sults  of  all  such  plots  made  by  tbs  SOS  are  summarized  in  Appendix  C. 

In  using  the  SC3  method,  c'urve  nu'bers  (Cu)  are  determined  by  re¬ 
ferring  to  a  table  such  as  that  sho'.?n  in  Appendix  C,  and  selecting  the 
CN  that  corresponds  to  the  appropriate  soil  type  and  vegetal  cover  con¬ 
dition,  The  soil  classification  used  divides  all  soils  into  four  groups 
(A,  B,  C,  and  D)  according  to  their  pemeabilities ,  Tne  permeability  of 
group  A  soils  is  highest  and  of  the  D  soils,  the  lowest.  The  descrip¬ 
tions  of  the  various  soil  types  are  given  in  Appendix  3,  Musgrave  and 
Roltan,  as  reported  by  Chow  (I96ir),  have  quantified  infiltration  rates 
for  the  four  soil  groups.  Minireum  infiltration  rates  for  groups  D,  C, 

B,  and  A  are  given  by  the  ranges  0  to  0.05^  O.Op  to  O.lp,  0,15  to  O.5O, 
and  0.30  to  0.45  inches  per  hour,  respectively. 

A  vegetal  cover  condition  Is  described  as  some  natural  or  culti¬ 
vated  conditicn.  For  instance  the  cover  my  be  described  various3.y  as 
fallo;;,  contoured  smiJL  grain  in  poor  condition,  pasture  in  fair  con¬ 
dition,  or  woods  in  good  condition. 

The  effect  of  the  anteccvient  noistiue  condition  (A.hC)  is  accounted 
for  by  adjusting  the  soil  cover  complex  number  previously  calculated, 

Tlie  curve  number  was  originally  developed  with  average  A.'iC  being  assumed. 
Kovever,  for  conditions  of  very  dry  and  very  v?et  antecedent  soil  moisture 

conditions,  type  I  and  type  III,  respectively,  are  used.  The  type  II 
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condition  is  considered  to  be  the  average  corditicn.  VJhether  the  AKC 
calls  for  type  1  or  'lype  III  treatreent  is  daterni>iel  by  the  guide  lines 
in  Table  1*  Tae  lir/its  on  the  dormant  season  apply  to  iinfroF.en  ground 
and  no  snow  covers 

After  ttie  cori'cct  type  of  AMC  is  determined,  there  are  standard 
curves  showi’ig  the  adjustments  to  be  Eiade  on  the  soil  cover  complex 
numbers  as  shown  in  Figure  2^  Valvies  for  the  adjusted  CN  can  bo  talxen 
at  points  batw’ecn  lines  shown  for  types  I,  II,  and  III  ':f  it  is  felt  the 
the  AFC  are  siiff icicntly  v’ell  hnown. 

St^riJ  Runoff;  The  sccon:!  subsystem  (that  which  produces  ths  storm 
runoff  function)  is  assumed  to  be  a  linear  convolution  by  unit  hydro- 
graph  theory*  The  purpose  of  the  cGnvo3-Ution  process  is  essentiaJl.y  to 
convert  a  ra.iijfall  excess  function  of  a  certain  volume  into  a  runoff 
function  having  the  same  voiAmiee  The  convoiution  integral  describing 
this  runoff  function  c](t)  is  given  by; 


q(t)  ::  r  ll(t  -t')  i  (r)  dr 
0 


(13) 


The  kernel  function  u(t  «  'b )  is  the  instantaneous  unit  hydrograph  (the 
unit  hydrograph  which  theoretically  results  from  a  ralruall  excess  occur 
ring  lnstantar'eou3.1y .  on  alJ.  parts  of  a  \7aterched ),  i{T)  is  the  rainfall 
excess  input  function,  and  t'  (the  limiting  time  f or  r )  is  a  function  of 
to,  to  being  the  duration  time  of  the  input  function  (storm)*  For 
t'  -  t  and  when  t  h.-  tQ,  t’  r  to* 

The  basic  ccr. volution  process  is  shown  graphical.ly  in  Figure  3* 
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Starbing  from  the  excess  function  i(t')  versus  ,  each  increricntal 
volume  of  rainfall  Is  multiplied  by  the  inrtantaneous  Uli  (HIH)  and  is 
suimned  in  the  right  time  sequence  to  yield  the  sto>m  runoff  function 
q(t)  versus  tc  The  instantaneous  unit  Jiydrcgraph  shoun  is  that  used  by 
the  SCS  when  developing  a  synthetic  imlt  hydrograph  for  an  irngaged  water 


ext^ 


shedc  This  Gynthetlc  UH  was  developed  for  the  SCS  by  Victor  lioclras 
reported  in  the  SCS  National  engineering  Eandbook  (196^1 )«  Olie  triangular 
lUH  sliown  \ritb  a  dotted  line  is  normally  uc-ed^  because  of  its  simplicity 
for  runoff  calculations  made  by  harle  However^  when  the  computer  is 
utilized,  the  mors  exact  ciu’v  ill  rear  lUH  is  eouai.lv  easy  to  annJr/e 

Trie  rainfall  excess  function  i(T)  can  be  used  directJ.y  with  the 
lUH  to  produce  the  runoff  function  q(t)<.  IIovever_^  the  v?ay  in  I’hich  val.ues 
of  i(C'  )  should  bo  calculated  is  not  obvious  from  the  earlier  develop- 
mont  of  equation  (3)«  Ihe  volume  of  rainfall  excess  i  (which  is  equi¬ 
valent  to  the  volume  of  runoff  q )  for  the  time  increment  dr  is  equal  to 
the  value  Q  given  by  equation  (3)  at  time  Z  dZ  minus  Q  by  equation  (3) 
at  time  Z  ®  The  volume  of  precipitation  P  is  the  total  precipitation  of 
the  storm  up  to  time  ft  cl?"  and  Z  ,  rcspectivelyc  Tiio  inoremenb  of 
(dQ);  caused  by  the  storm  of  df  length  must  be  given  as  a  difference  of 
Q's  at  times  ft  d  f  and  Tais  is  true  since  the  runoff  equation  is 

based  on  the  entire  storm  (the  inibial  abstraction  must  be  accounted  for 
only  once). 

The  base  time  for  the  curvilinear  UII  in  Figiu’e  3  is  taken  to  be 
times  the  time  to  peak.  This  value  of  base  time  ifi  arbitrary  but  docs 
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nci;  incorporate  appreciable  erior  in  tha  raoirlts  as  the  ccrresponvling 
values  of  discharge  are  very  sl-sII,  Ihe  peak  flow  rate  for  the  cur“ 
vilinaar  UH  is  the  same  as  that  for  the  triangular  vex s ion  uhcss  base 
tii/ie  (determined  mathciratlcally )  is  found  to  be  2,67  times  the  time  to 
peak  (tp),  Ihs  triangular  liydrogi’aTjh  base  time  of  2.67  ip  is  used  as 
it  gives  the  same  value  of  runoff  voluine  as  the  cui'vilxncar  version^  the 
volume  being  repreGented  by  the  area  under  the  curve.  The  iDeak  flow 
rate  Op  is  derived  using  the  triangular  shaped  UII  and  is  given  by  the 
ge  no  ra  1  e  q  u.at  i  on 


ip 


'P 


\?ith 


tp  s  d/2  L  K  I)/2  1  0,6  i'Q 


(15) 


a  no. 


tg  “  2,67  tp  for  triangular  unit  hydrograx/ns 


(16) 


and 


t’jj  r  5*0  t  for  cm’vilinear  unit  hydregraph: 


(17) 


uhere  K  Is  a  constant  whose  value  dependG  on  the  mrits  used^  A  is  the 
arca^  Q  is  total  volume  of  runoff^  tp  is  time  to  pss.k^  is  the  base 
time,  I)  is  BtorBi  duration  time^  tc  is  the  time  of  concentration^  and  I, 
is  the  basin  lag  tl-vei  K  is  equal  to  kQh  when  0  is  in  incheS;  A  in 


SQuarc  miles  and  t,,  L  end  D  are  in  houins 

P' 


Peak  flow  is  then  in  cubic 


feet  per  second 


For  the  example  shown  in  Figure  3;  where  the  convoluticri  integral 
gencsrally  cannot  be  solved  exactly^  actual  px-actice  dictates'  that  the 
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storm  Input  ba  broken  up  into  srfiull  Inerenants  a  r  .  Tne  storiii  duration 
D  of  eq\iaticn  (I5)  is  then  equal  to  a  r  *  Equation  (l4)  toscther  vrj  th 
either  (3„6)  or  (17)  is  sufficient  to  describe  the  rimoff  hydrograph  from 
a  subuatershed  due  to  an  incremental  storra.  The  liydrograph  can  be  given 
the  triangular  or  curvilirear  shape  as  desired. 

F.lood  Runoff ;  The  third  subs^/stem  (that  uhicb  produces  the  flood 
runoff  function)  is  the  simple  procedure  of  adding  the  baseflow  hydro- 
graph  to  obtain  the  resulting  flood  hydrograph.  T.as  baseflou  hydregraph 
may  be  obtained  by  any  of  the  standard  textbook  methods  as  long  as  the 
same  procedure  is  used  consistently. 

Channel  Floci  Rout ing 

V/hen  a  watershed  is  divided  into  subvatersheds  for  the  purpose  of 
predicting  runoff  more  accurately^  it  is  necessary  to  route  hydi^ograxfns 
through  the  natural  channels.  Tlie  routed  bydregrax^h  (outflow  hydro¬ 
graph  from  a  channel)  is  then  in  a  proper  form  to  be  added  to  another 
routed  hydrograph  or  to  another  subwaterslied  runoff  hydrograph,  after 
which  it  may  be  routed  through  a  second  reach  of  channel, 

Fnc  method  chosen  by  the  SCS  for  channel  flood  routing  is  termed 
the  "convex  method"  after  the  name  applied  to  the  mathematical  principJ.e 
upon  w'hich  it  is  based. 

The  convex  method  is  based  on  tVie  principle  that  when  a  flood  passes 
through  a  reach  of  natural  channel  there  is  a  time  interval  such  tliat 


if  Ij  kb  Op,  then  Ij  .V  0^  kb  0^^ 


(V'a) 
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if  0^,  then  02  0-^ 


(It 


) 


yhcrc^  J.j_  aM  0^  are  'cho  rates  of  inflo:.’  and  outfloy  respectively  et  time 
t  ::  t2,  and  Op  is  the  rate  cf  outflow  at  time  t  r  s  tp*-  Note 

l-hat  Dp  (D.uflow  rate  at  tp )  does  not  enter  into  this  method*  If  At  is 
chosen  as  specified  then  Oj ,  and  Op  form  a  "convex  set".  From  con¬ 
vex  set  theory  we  can  tlien  write  the  uor-rkinj  equation 


Op  (1  -  C)  0^-bCI^ 


(IS) 


v.’herc*  C  is  a  parameter  with  the  range  O-C-'l. 

From  the  definition  of  a  convex  set,  it  can  be  said  that  yhanover 
any  two  points  x  and  y  are  in  a  convex  set  L^the?!  every  point  on  the 
J.ine  connecting  them  is  also  in  L»  A  formal  definition  from  set  theory 
given  by  Charnas  and  Cooper  (I96I)  states  that:  a  set  of  points  L  is 
convex  if,  w’ he  never  x,  y-^L-,-  then  u>l  -h  (J.  -  u)y-sL  for  all  O-u-le 

The  Units  on  the  value  of  C  insure  that  equation  (19 )  meets  the 
conaitions  of  the  equations  (l8a)  and  'l8.b).  It  follows  that: 


C 


tl 


0.1 


(20) 


It  can  also  be  sho'ni  frora  Figure  t  that: 


O2 


0 


1 


-■Cl  “  0l 


A  t 


(21 ) 


v;hcrc  K  is  a  parameter  '.?ith  time  units  not  yet  having  any  physicaD  mean* 
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Figure  4:  Fg  la  t  ions  hip  hotv/een  At  and  K 


ing. 


From  equa-Giori3  (20)  and  (21)  it  can  be  seen  that 


h  t 


CK 


II  use  of  this  method  is  to  bs  rade  on  ungaged  o’atsrsliods ,  soma 


empirical  re lat lonships 
K.  This  has  bean  dons 


need  to  be  developed  for  evaluation  of  C  and 
by  the  SC3  from  knowledge  gained  from  study  of 


gaged  uatorsbeds. 


It  was  found  that 


equation  (23 )  expresses  the  value 


C 


(23) 


for  C  very  veil  except  at  extreme  values  of  V(<0,h  and  >-10*0) 


where 


V  is  the  average  stream  velocity  In  feet  per  second 
bo  the  travel  time  for  the  roach.  If  K  is  in  hours 
le rigt h  in  f  e a t  ^  the n  i 


,  K  \?as  also  found 
and  L  is  the  reach 


to 


K 


3600  V 


(.2>: ) 


It  is  obvious  from  the  empirical  relationships  for  C  and  K  that  as 
the  v?atcr  surface  c3.evatlon  changes^  varying  velocities  will  cause  C  and 
K  to  vax-y  vrith  tine.  A  procediu’c  for  determining  V  W’as  sought  so  that 
C  and  K  could  be  assumed  constant  for  any  given  routing  ])roblcm.  It  was 
found  that  if  the  average  vel-ocity  V  was  based  on  the  stream  flow  velo¬ 
cities  for  all  discharges  of  the  inflow  hydrograph  greater  than  ono-kualf 
of  the  peak  discharge;  the  resulting  outflow  bydrograpb  was  nearly  the 


Dame  as  if  it  were  calculated  using  varying  C  and.  K, 


i 
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Tiae  va3Aio  of  At  to  ba  used  is  then  a 
uiains  constant  for  the  entire  routing,  The 
to  be  iDcoaveiiient  fractions  of  an  hour  and 
SCS  practice  is  to  choose  a  convenient  ' 
C',  such  that  routing  by  equation  (19)  give 


product  of  C  and  K  and  re¬ 


value  of  At  nay  turn  out 


the  re  f  ore  d  iff  ic  u].t  t  o  us  e  « 
and  calculate  a  neu  C.  called 
s  identical  results.  The 


acijustm-ent  is 


1  ^  (1 


j  At  '  /  At 


(25) 


Routing  is  nov/  accomplished  by  rewriting  equation  (I9)  to  give  the  result 


Op  »  (1  “  C  )  Oj 


Reservoi2" 

Runoff  h^/drographs  are  routed  through  reservoirs  by 
dication  method.  This  n.ethod  is  based  on  the  continuity 
here  in  its  differential  forn_j 


the  storage  in¬ 
equation  sho'.’n 


i(t ) 


(^87) 


where  i(t)  is  the  rate  of  the  inflow  at  time  t,  q(t)  is  the  lafce  of  out- 
flow  at  time  t,  and  dS/dt  is  the  rate  of  change  in  storage.  The  basic 


v;  or  king  equation  is  written: 


(Il-Mp)  At/2  (0^-{-  0^)  At/2 


Gp  -  Sj_ 


(28) 


w’here  I  is 
subscrints 


inflow,  0  is  outflow,  S  is 
1  and  2  are  for  t  t-j  and 


storage,  t  is  a  time  increment, 
t  s  t  j_  +  r.t  s  tp  respect  ivoJ.y, 
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Keduction  of  this  gives; 

( -r  I2  "!'■  K  Op  .4,'b  V  2 So 

DafiDing  constants  and  as  : 

C4  ~  ( I  -{-  Ig  -  0i )  t  -1-  2S^ 
02i^i‘t'r2S2 

It  then  fcdlons  that  fi’o-a  eqijation  (29)  J 

C4  =  C5 


(29) 


(30) 

(31) 


(32) 


C]^  is  first  evaluated  from  knoun  values,  then  set  equal  to 
This  value  cf  C-  is  then  substituted  into  equation  (31)*  There  is  one 
set  of  values  (0^  and  S2  )  on  the  storage  “discharge  cuTve  for  the  reser** 
voir  that  satisfies  this  linear  equation.  This  procedure  constitutes 
.one  step  in  solving  the  continuity  equation  by  nuraerlcal  approximation. 
To  route  a  bydrograph  through  a  reservoir  this  process  mu.st  be  repeated 
until  flood  flous  have  nearly  dininished  to  the  original  basoflov/. 
Systems  Class  if  Icat  ion  o_f  S^S  Method 

The  runoff  dynaraics  cf  an  actual  vatershed  night  be  described  as  a 
distributed,  non-linear,  time  variant  system,  Ihny  of  the  models  con¬ 
structed  to  describe  this  system,  however,  are  c-lassificd  as  lumped, 
linear,  and  time  invariant.  Early  models  \;ith  these  assumptions  v?cre 
those  devc.loped  by  Sherman  (1932)  and  Bernard  (J.935)^  Sherman  developed 
the  unit -graph  concepts  often  called  unit  bydrograph  thaory.  Bernard  is 
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given  the  credit  fcr  iiupJ.enient ing  use  of  the  distr5-hut ion  graph  (a  bar 
graph  comparable  to  a  unit  hydrograph  cncept  that  runoff  volumes  are 
expressed  as  percent  of  total  runoff  versus  tine), 

A  distributed  system,  in  the  field  of  hydrology,  is  defined  as  one 
with  time  and  space  variation  in  rainfall  intensities,  A  lumped  system 
or  model  in  one  which  allows  only  time  variation  but  no  space  variation 
in  rainfall  intensity,  Ihis  is  to  say  that  rainfall  intensity  is  uniform 
over  tJi.e  watershed  area  under  consideration  at  any  given  time. 

The  SCS  method  in  Its  basic  forra  (as  currently  explained)  would  be 
classified  as  a  luiTipcd  systemc  By  using  the  SCS  i-unoff  prediction  method 
in  conjunction  with  resei’voir  and  channel  flood  routing  techniques,  the 
resulting  system  ^’i  11  in  part  be  distributed,  Chow  (1967)  describes  this 
type  of  system  as  a  "distributed  system  of  lumped-'system  models".  On 
each  of  several  subwater  sheds  the  rainfall  intensity  is  uniforra,  or  each 
of  the  subsystems  is  lumped.  However,  at  any  given  time,  the  intensity 
can  vary  from  subw-atershed  to  subwatershed  and  is  in  this  sense  distri- 
buted. 


A  system  is  linear  if  the  unit  responses  of  various  input  iragnituies 
are  identical.  The  watershed  runoff  function  is  linear  if  one  unit  hydro- 
graph  is  sufficient  to  accurately  describe  all  storms  of  a  given  dui’ation 
occurring  on  a  given  watershed.  This  is  never  the  case,  however,  as 
Childs  (I95S)  and  others  have  shown.  Childs  states  tViat  the  larger  the 
flood,  the  higher  the  peak  and  the  shorter  the  time  to  peak  on  the  unit 
hydrograpb.  He  demonstrates  this  fact  by  showing  four  unit  hydrographs 
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of  three -hour  c-torrr.  on  Faugatuck  River  at  Thoirastori;  Conruectic^it . 
The  peak  of  a  ina jor  flood  baa  raore  than  double  the  peak  derived  for 


a 


minor  flood  and  about  tuo-tlri.rds  of  the  tin^  to  peak,  Cbvloualy;  the 
variance  of  the  storm  pattern  or  rainfall  excess  function  af'fects  the 
unit  hydrograph  shape.  ITon-llnearity  is  then  the  case  uhen  the  derived 
unit  hydrograph  varies  uith  storm  pattern  and  total  rainfall  excess  amount 

A  system  that  gives  the  same  response  for  a  given  input  regardless 
of  time  is  a  time  invariant,  system.  The  subuat-ershsd  portion  of  the  SC3 
model  is  based  on  unit  hydrograjsh  theory  and  is  therefore  tira-e  invariant. 
This  is  not  too  restrictive  an  assumpirlon  as  the  rairl'all  input  to  the 
Eubu’atershGd  riinoff  moclel  is  preadjusted  for  soil  and  cover  conditiers; 
antecedent  iaolsture,  and  time  of  yearo  These  are  the  main  items  giving 
tlie  runoff  fu.net ion  (system  response)  a  time  variant  iintura. 

Spatial  and  temporal  variations  in  rainfall  distribution  can  be 
accounted  for  vritb.  the  SCS  apxjrcach.  These  variations  are  simulated  by 
dividing  a  uatersbed  into  several  suboatersheds  and  by  making  the  time 
increment  on  the  precipitation  versus  tine  input-  gi’a^fri  small  enough  so 
that  the  gx^aph  nearly  matchSvS  the  actual  curved  relationship.  The  total 
amounli  of  rainfall  can  be  different  for  each  subuatershedc  The  storm 
pattern  can  have  an  infinite  number  of  shape  variations.  Spatial  var- 
ration  is  not  totally  accounted  for  since  there  is  still  no  spatial  var¬ 
iation  allov?ad  uithln  the  subwatorshsd ;  thus  spatial  variation  is  only 
partially  simulated.  The  cubw'atersheds  can,  hov?ever,  be  made  as  small 
as  is  practical  for  the  time  and  expense  that  can  be  spent  in  the  ara^-V-sis 
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T/ie  tlina  of  storra  boginnin.?:  and  duration  can 


necessary  to  confoinA  to  an  actual,  storin*  Therefore 


Iso  be  designated  as 
tenporal  variation 


is  also  partla.lly  siraulated^ 
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Chai>ter  lY 


SOIL  COESEHVATIOII  SI-RVICE  PROJECT  FOKidULYTiON  PROGRAM  -  HYDROIOGY 
The  Soil  Confiervaticn  Sarvica,  in  1963;  cbtained  (froLi  a  private 
conouj.tlri^  fiXid  kno'.-n  as  R*,  Ins.)  a  digital  computer  program 

designed  to  solve  tlie  ’’runoff  hydrograph  prediction''  relationships  dec— 
cribed  in  Chapter  III*  The  computer  prograL'U  v/hlch  uas  vrittan  in  the 
Fox’tran  language  for  the  _  I}3M  709O  computer  system^  is  entitled  'lirojcct 
Formulation  Program  -  Ilydrolcgy"  and  is  a  rathairuitical  reodel  of  a  uater- 
shed  experiencing  precipitation;  surface  runoff;  and  channel  flc.?.  The 
program  can  route  an  unlimited  number  of  hypothetical  or  actual  storms 
through  a  watershed  having  as  many  as  60  structures  and  IfO  crossoscc- 
tions  described o 


0?he  program  has  gone  through  several  revisions;  most  recently  in 
October;  1967 >  It  is  non  uritten  in  Fortran  IV^G  and  utilises  the  IBM 
360/40  comiuiter  system.  The  most  recent  instructions  for  using  the 
PTvOject  Formulation  Program  (PFP)  uerc  published  by  C.  E,  I.  R.;  Inc, 
(1964)  and  the  SCS  (1965), 

T3]e  PFP  consists  essentlaJ.ly  of  a  series  of  subroutines  uhicli  ray 
be  so.lved  in  a  variety  of  combinations  and  sequences;  a  schematic  over¬ 
view  of  the  PFP  is  shcrni  in  Figure  5*  Thu  particu-lar  sequenca  to  be 
folRou'ed  in  a  given  problem  is  detcriuj-ned  by  instructions  to  a  "Stan¬ 
dard  Control"  routine;  vdrich  calls  the  subroutines  in  the  Xjroper  order; 
and  provides  the  additional  information  neces.sary  to  perform  the  work 
of  the  subroutine.  The  entire  program  is  juonitored  by  an  "Executive 
Control";  v/bich  describes  the  x-^i’^icnAlar  stcrin  under  coiiDideration, 
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To  fonnulat-e  a  ?PP  problem;,  it  ia  first  zieccssary  to  read  in 
several  tables,  as  ind-Jcated  on  Figure  5^  soro  general  and  somu  spa- 
cific  for  the  particular  watershed,  -  After  this  is  dons,  Standard  Con¬ 
trol  cards  are  a’oad  in  to  describe  the  logical  sequence  in  which  flood 
routings  through  tiie  stream  reaches. and  structures  should  be  performsd, 
•This  sequence  of  opaiations  is,  cf  course,  identical  to  the  sequence 
that  the  hydrologist’  follows  in  solving  the  problem  cianualljc 


To  reflect  a  future  change  in  a  watershed,  ary  portion  of  the 
Standard  Control  can  be  modified  at  any  time  by  addition,  alteration,  or 
de.letion  cf  \jatersli3d  data.  Executive  Controj  cards  are  than  read  'Jnto 
the  coiaputer  to  dascriba  each  of  the  alternative  storm  situations  that 
are  to  be  analysed.  An  infinite  number  of  problems  may  be  set  up  and  run 
by  the  Ebiecutive  Control  for  a  given  set  of  watershed  characteristics. 

Each  directive  of  the  Executive  Control  specifier  the  storm  precipi¬ 
tation,  its  starting  tiir.a,  th.c  antecedent  moisture  condition  under 
which  it  is  to  be  analyzed,  and  the  portion  of  the  watershed  through  which 
it  is  to  be  routed, 

Ka.ch  sub'/etershed  described  may  have  different  total  rainfall  amounts, 
different  storm  starting  tinics,  difforenl:-  storm  tine  dm-ations  and  pat¬ 
terns,  and  varying  antecedent  moistui’e  corditiou.s,  A  given  storm  sit¬ 
uation  is  then  applicrl  to  the  Standard  Control  to  effect  a  so.lution. 


Additional  storms  may  then  be  set  up  and  solved.  As  has  been  seen,  the 
Standard  Control  can  at  any  time  he  changed  to  reflect  propose;!  changes 
on  the  watershed.  After  a  charge,  any  hypothetical  or  actual  storm  can 


again  be  applied  to  the  modified  watershed  to  evaluate  the  chanses  made 
3.n  the  outflow  hj/drograph. 


■  The  SCS  program  has  a  total  of  twelve  GubroutinsG  (seve’n  of  which 
are  computational ) o  Their  relationships  to  the  inai^n  program  ani  to 
other  subroutines  are  shown  in  Figure  6o 

The  seven  corrout  at  tonal  subroutines  are  called  EICOFFX,  R.S3V0Ii_, 

REACH;,  ADDini)^  OUIPlFi?^  IHTERP,  and  PEAKc  They,  respectively,  calculate 
a  runoff  hydrograph  for  a  cubwatersbadj  route  an  input  hydrograph  through 
a  reservoir  by  the  st-orage- indication  method;  route  an  Input  hydrograph 
through  a  channel  reach  by  the  convex  method;  add  two  hydrogruTdis;  print 
and/or  punch  the  output  options  specified;  adjust  the  time  increment  for 
RESVOR,  REACH,  and  ADDilYD;  and  calculate  coordinates  of  the  highest  peals 
(to  a  maximum  of  ten)  of  a  hydrography 

The  other  five  subrout inas  are  called  SAVMOY,  UPPATE,  SETUP,  READ  5, 
and  READINo  SAVMOV  gives  the  v;riter  of  Standard  Control  the  ability  to 
move  an  entire  hvclrograph  from  one  storage  location  to  another  and  to 
save  it  there  until  needed  later*.  UPDATE  causes  the  Standard  Control 
listing  currently  in  force  to  be  printed  out*  This  would  be  used  after 
some  alteration  of  Standard  Control  had  bee'u  made*  SETUP  has  the  over¬ 
all  control  of  the  coraputational  subroutines*  PiiA.D  5  RHIADIN  cause 
the  program  to  read  in  the  data  \.’hich  describe  the  physical  features  of 
the  v/atershed  and  the  properties  of  the  storm  being  simulated* 

Tabular  I'ata 

Several  t>pes  of  data  that  are  necessary  for  the  operatior  of  the 
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Project, 


SclieiTie  for 
ForrAilc.t  Ion 


Soil  Conservation  Service 

i  X  V  (-W-L  V..!. 
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Project  Forniulation  Program  can  be  Bup^^lled  to  the  coupirtsr  in  tabular 


3?orm.  These  data  are  valnss  for  the  relationship  botweon  the  routirig 
coefficient  C  and  average  ve.locity  V,  the  curvilinear  n.nit  hydrography 
sterm  patterns,  the  stage-storage-clis charge  relation  for  reservoirs  and 
the  stage -area  “discliargo  relation  for  channel  crosS“Section.' ^ 

The  C  ”  V  ga'aph  is  tabulated  for  co:aputer  use  in  velocity  incre- 


in  t/tp«  resign  storiTi  patterns  aie  described 


rronts  of  0,2  feet  per  seconds  The  unit  hydrograph  is  supplied  in  dinen 
sionless  form  as  cj/q 
in  cumulative,  dii iansicnlcss  form  uith  coordinates  of  P/p.p.  and  t/tp 
vhere  Fp  is  the  Gtorm's  total  procip>itatlon.  Actual  storns  can  be 
read  in  either  dimensioned  cr  dimension.lc ss  form  but  sluays  as  a  cumu¬ 
lative  graph.  If  actual  dlr.ensions  are  used,  rainfall  depth  and  duration 
in  Executive  Control  are  to  be  given  as  1,0.  If,  heuever,  diiaa  ns  ionless 
coordinates  are  used  for  actual  stornis  in  tabular  data,  the  actual  values 
of  depth  and  dirration  irast  be  supplied  in  the  Executive  Control.  The 
structure  and  cross-section  data  can  be  given  in  any  dcsiiablo  increr.ient 
of  elevation,  Ihls  should  be  a  suall  enough  increment  so  that  the  assured 
straight  line  betuoen  points  is  nearly  coincident  yj.th  the  ciu’ved  re¬ 
lationship  from  vfnich  the  points  v^ere  selected,  The  elevation  incre¬ 
ment  can  be  varied  at  uill  throughout  these  tables. 

The  vay  in  uhich  these  t8.bular  data  are  utilised  by  the  Executive 
and  Standard  Control  Is  outlined  by  the  flov?  chart  in  Figure  5^  lire 
V)rob.lems  posed  by  the  Executive  Control  are  individual,  distinct  storms, 

■an  infinite  nurfocr  of  uliich  can  be  processed  by  the  PF?«  After  probiem 
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Ko,  1  hag  been  vjorked  by  tba  Btandard  CoBurol  through  jt-s  various  r.ub- 
routines,  problem  NOc  2  ic  then  described,  by  the  Executive  Control  aril 
solved*  This  continues  until  aJl  stated  problems  Irave  been  solved* 

These  subroutines  in  turn  use  the  tabular  data  (shown  in  squares  on 
Figure  5)  complete  their  function*  REACH  needs  the  C  -  V  and  the 
cross-section  graphs  to  complete  its  task.  The  cross-'sect ion  information 
is  needed  whether  Ou’TPUT  is  called  or  not.  RI'IOl'R’X  uses  tb.e  UH  and  either 
a  design  or  actual  Btorm  P  -  t  gra'ph.  SAVIiOV  and  ADDHiD  do  not  make  use 
of  the  tabular  data*  RS3V0R  must  calJ.  the  structure  data  to  perform  its 
routing  problem.  All  five  of  the  working  subroutincfi  may  hc.ve  to  call 
OUTPUT  (if  specified  by  standaixl  control)  before  returning  to  Standard 
Control. 

Standard  Ccm^rol 

Standard  Control  is  an  ordered  set  of  instructions  to  the  computer 
that  outline  the  sequence  of  operations  that  simulate  runoff  from  a  water¬ 
shed  area.  Tliese  operations,  which  have  aJ.ready  been  outlined,  are  done 
by  the  subroiPrinsG  RNOPFX,  RE3V0P,  RF^iCH,  AL'DIIYD,  and  SAVMOy.  PU.gure  7 
sho-.-'S  a  sample  watcr?Rn&d,  with  two  reservoirs  and  two  main  stream  branches, 
which  is  subdivided  into  five  subwaterslieds .  /\  tabi.e  of  the  correspond- 

irig  Standard  Control  listing  for  this  sample  vmtor-shed  is  also  given, 

Ihe  structure  and  cross-section  locations  are  nuribered  so  that  the  reader 
can  follow  through  the  Standard  Control  to  see  that  iu  adequately  describes 
the  v;ay  in  wliich  rainfall  excess  is  converted  into  runoff.  There  arc, 
however,  soma  assumxjticns  in  this  connection  wliich  make  the  sirimlaticn 
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Figure  7  i-  SarupJ-C  V/aterohed  u5.th  Correspoiiriing  Standard  Gonti'Ol  for  the 
Soi.l  Conservation  Service  Project  Forinujat ion  Program 


ihMch  easier  ulthcui;  acui5j:g  substantial  error  to  the  operation.  For 
exaiaplG^  the  out-flou  hydrograph  from  structure  1  is  roubed  tbro',-\gh 
area  II  ao  if  no  addit5.onal  runoff  v’ere  being  supplied  from  area  H, 

Nov_,  there  exists  an  outflow’  hydrograph  froHi  the  channel  reach  in  area 
II  at-  section  1,  To  this  hydrograph  is  then  added  the  runoff  liydrograph 
from  area  II  as  if  it  had  not  arrived  at  section  1  by  way  of  the  channel 
system.  This  assumption  is,  of  course,  not  physically  the  case  but  soon 
after  Dmaving  the  upstream  subwatershods  the  channel  hydrograph  is  very 
large  compared  to  the  local  area  runoff  so  that  the  error  resulting  from 
this  assumption  is  considered  scmll. 

Computer  locations  for  hydrogra’ph  storage  are  shown  on  the  table 
included  in  Figure  7c  There  arc  seven  different  locations  sux^plled  where 
all  the  coordinates  of  points  from  a  bydrogranh  can  be  stored  (up  to  300 
points),  Tiiese  x^oints  are  all  described  by  a  time  and  dis^charge  value. 
All  subi^outinss  have  an  input  and  an  output  hyd.rograph  except  ADDIiYD  vdiicl 
has  two  input  an-d  one  output  hydrograjphs ,  The  numbering  of  storage  lo^ 
cations  as  shown  is  set  and  consistently  used  exce]pt  for  the  SAVMOV  sub¬ 
routine  which  can  then  be  used  to  place  a  hyd.rograph  in  the  correct  loca¬ 
tion  for  the  subsequent  operationo  This  rigid  numbering  system  is  used 
60  that  fewer  errors  will  be  made  by  a  technician  preparing  the  Standa.rd 
Control. 

All  the  giossiblo  output  options  are  also  shown  in  Figiu’e  7*  These 


include  printing  out  all  the  values  of  the  hydrography  the  corresponding 
water  surface  elevations,  up  to  the  ten  highest  pc-nk  discharge  values  and 


their  timec  of  oecurrtjncc*,  and  the  total  volru-.e  under  the  bydrograpb. 
given  in  three  sets  of  units,  cfs»hrs,  acre^feet,  and  inches  of  runoff 
over  the  watershed <> 

Obher  informtion  than  that  shown  must  also  hs  supiJliecu  Area 
in  square  miles,  runoff  cervo  nuriher  CN,  and  tire  of  concentration  in 
hours  riust  aJ^vOO  be  given  with  the  KiJOFFX  cammar-d .  HESVOR  must  also  knoi? 
the  surface  e.lGvation  of  the  I'eservoir  at  tim.e  t  s  0.  length  of  a  roach 
In  feet  and  a  routing  coefficient  (if  one  is  to  be  supplied  instead  of 
being  calculated  by  the  program)  are  needed  to  complete  the  P.EACIi  Gub-» 
rou.tinec  ADDliYI)  and  SAVMOV  are  complete  as  shown  in  Figure  7® 

Viith  regard  to  the  numbering  of  stream  channel  cross -sect ions  and 
control  struct ui'es,  numerical  values  my  be  given  in  any  order,  a.lthough 
a  logical  sequantial  order  is  easier  for  others  to  follow  *  However,  the 
\v’ay  in  which  these  numbers  are  used  in  Standard  Contz’ol  is  very  important 
as  will  be  seen  when  the  use  of  the  Executive  Control  is  explained.  Tne 
important  thing  to  remember  Is  that  any  number  used  more  tho.n  once  must 
appear  consecutively  in  the  Standard  Control  listings  That  is,  if  cross*’ 
section  001  is  used  as  it  is  in  Figure  7^  four  times  that  it  appears 
must  follow  one  another  without  interruption  by  other  nvimhero. 

Executive  Control 

Executive  Control  is  a  set  of  comriinnds  •  that  liave  the  over  “all  con¬ 
trol  of  the  coraputcr  solutioji  to  stated  hydrologic  probixzns.  As  in¬ 
dicated  earlier,  each  Executive  Control  directive  specifics  the  proper* 

tics  of  the  iiypothet leal  or  actual  storm  that  is  to  be  routed.  Storm 
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propertiec  accounted  for  in  this  Riodel  Include  &torm  pattern^  total  rain¬ 
fall  amount,  otorm  length,  and  starting  time.  Tae  proparties  r.ay  vary 
frora  ouowatershed  to  subuaterched  but  are  assumsd  uni.forn  over  each  in- 
d  iv  idua  1  s  ubwat  e  rahe  d  <, 

The  storm  pattern  is  generally  read  into  the  computer  with  a  dimen- 
Eionless  cumulative  raiiifall  amount  ordinate  and  a  dimensioned  time  co- 
ordinate c  However,  it  is  pcssihle  to  read  in  both  coordinates  with  either 
dimene-ionsd  or  dimensionl.ess  scales.  Tne  dimensionless  values  are  then 
transformed  to  dimsnoioned  values  by  specifying  the  total  rainfall  amount, 
the  total  storm  length  that  applies  to  the  particular  case,  or  both.  This 
flexible  arrangement  is  useful  for  applying  hypothetical  storms  to  the 
system.  In  the  case  of  actual  storms,  cumulative  graphs  with  tlr'S  correct 
time  and  precipitation  diffiensionG  can  be  read  in  directly. 

INCREM,  COlfPUf,  ENDCiiP  are  additional  statements  used  in  Executive 
Control  which  simply  are  shorthand  for  increment,  compute,  and  end  of 
computations,  respectively.  They  are  not  subroutines  but  simply  commiands 
that  are  used  by  the  main  program  to  carry  out  the  desired  control  se¬ 
quences  « 


Examples  of  Execiit ive  Control  t  Several  example  xxroblems  are  iJ,- 
lustratecl  in  Table  2  to  dofflonstrate  the  flexibility  of  the  Executive  Con' 
trol.  Problem  Eo,  1  breaKs  the  v/atershed  into  three  subwatersheds  each 
receiving  different  amounfcs  of  rainfall.  Rainfall  Treble  No.  1  would  be 
given  in  dimensionless  depth  units  and  time  in  hours  since  the  amount  of 
rainfall  is  specified  and  1,00  is  given  under  rainfall  duration.  The 
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Table  2:  Exaiiple  of  ExecutiYQ  Control  cn  Sample  Watershed. 


anbecedcnt  moisture  condition  (AMC)  is  type  II  or  normal*  Gtarbing  timo 
is  given  as  zero  on  all  su>jv;ater sheds* 

Pr-oblem  Ho.  2  is  similar  to  No,'  1  except  that  storm  movement  is 
accounted  for  by  varying  the  starting  tine  021  the  various  subv/aterehsds * 
Also  the  storm  pattern  is  considered  to  be  different  in  shape  or  lengbh 
as  rainfaj.!  Ifeble  No*  3  is  used*  The  AMC  is  type  •  I  or  very  ary ^ 

Problem  No*  3  is  an  actual  storra  that  is  taken  as  uniform  In  amount 
over  the  entire  v’atershed*  Its  starting  time  varies  for  tu’o  subareas  and 
a  nev.'  Tab3.e  (Ko*  2)  describes  the  storm  patterns  uhich  is  given  in  true 
units  of  inches  and  hours  since  the  factors  used  for  rainfall  depth  and 
duration  are  both  given  as  unity.  Problems  4  and  5  s2.-e  design  storms  and 
are  uniform  over  the  entire  area  in  rainfall  amount  and  starting  tine, 
llie  rainfall  dcxybhs  given  uouJ.d  correspond  to  some  calculated  design  storm 
amount  for  a  given  frequency  storm  for  the  area  in  question,  Painfall 
Tables  and  5  I'iight  be  storm  imitterns  for  12  houjrs  and  24  hour  storms. 


A  normal  AMI  \?oiild  always  be  assumed  for  design  storm 


INCREM  states  the  time  increment  in  hou.rs  that  the  ecmiputer  should 
use  in  developl2'g  the  hydrogra'phs ,  COI'T^uT  instructs  the  computer  to 
compute  what  is  called  for  in  vStanuard  Control  between  the  structure  (s) 
and/or  crc)ss«-scction(s )  indicated  in  the  command,  ENDCM?  indicates  the 
end  of  a  COhlPUT  series  and  sends  control  back  to  the  beginning  of  Standard 
Contro’L  to  start  a  new  runoff  problem_j  based  on  another  storra  described 


in  the  Executive  Control, 


ChaiTte?’  V 


PRiilPAimTIOw  0?  V/ATERSLED 


The  preparation  of  uatershed  and  storn  data  for  use  uith  the  5ioll 
Conservation  Service  Project  Formulation  Program  (PFP)  vrill  now  be 
demonstrated o  The  examp.le  used  vill  be  the  Duck  Creek  Watershed  near 
Brockway^  r-vontana  and  the  storm  of  jure  iG^  I965,  Figure  8  is  a  map 
of  the  Duck  Creek  clrainap;e  shoeing  the  location  of  rain  gages  and  rater  ' 
stags  recorders  o  The  stream  network  is  aI.EO  clearly  seen;  as  is  the  sub¬ 
division.  of  the  watershed  into  sub'satersheds . 

The  largest  of  the  four  watersheds  mider  study  by  the  Drainage 
Correlation  Research  ibroject^  Duck  Creek  (5^  square  miles),  located  south 
of  Brockv/ay,  Montana,  \ias  selected  for  the  investigation  reported  herein^ 
it  v/as  chosen  for  several  reasons:  l)  its  size;  2)  the  occurrence  in 
1965  of  a  rain-caused  runoff  event  v.’hieh  apxieared  to  be  suitable  for 
analysis;  and  3)  particular  interest  in  this  watarshed  by  the  Montana 
State  office  of  the  SC3,  Instruments  instaJ.lcd  at  Duck  Creek  by  the 
Drainage  Correlation  Research  Project  include  three  V/cather  stations 
u'hich  measure  wind  velocity,  air  and  soil  temperati.u’o ,  and  soil  moisture; 
four  non- recording  and  three  continuous -recording  rain  gages;  and  tv^o 
continuous  “recording  water  3-evel  recorders.  One  water  .level  recorder  is 
located  on  the  bast  Fork  of  Duck  Creek  where  it  is  crossed  by  State  High- 
v;ay  253 <>  Ihe  other  is  located  on  IXick  Creek  330  feet  north  of  the  center" 
line  of  the  county  road  between  sections.  11  and  it,  TI7H,  R46s. 

The  PF?  is  a  generalized  T.irograin  designed  to  be  used  on  a  wide 
variety  of  watersheds,  and  is  capab.le  of  processing  many  types  of  data. 
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AnaJ.ysis  of  the  Duck  Creek  v?atershecl  did  not  require  utilization  of  all 
the  program  capabilities  *  liie  data  supx^-jed  to  the  coiuputer  for  Duck 
Creek  relate  to  l)  cubuatershed  area's  and  channel  lengths^  2.)  stage- 
arcca  and  stage  “discharge  relationships,  3)  curve  numbers  and  antecedent 
moisture  conditions,  4)  time  of  concentration,  and  5)  Btcrm  data.  These 
are  discussed  in  detail  balou’, 

Subwatershed  Areas  and  Channel  Lengths  * 


For  'purposes  of  construct i.ng  a  model  of  the  Duck  Creek  channel  net- 


vork,  the  Hack  Creek  watershed  v?as  subdivided  into  twenty-four  au.bwat 


e:< 


sheds 


the  projected  area,  channel  meander  length,  and  the  flood 


length  Vrciwi  measirred  for  eacli  of  the  suiowatershads .  Areas  were  deternined 
by  p3,anircetcring  aerial  'photographs  (scale:  4"  u  1  mi,)  supplied  by  the 
SC3,  Both  insanider  and.  flood  plain  lengths  were  measured  on  the  same 
photograi)hs  v?ith  a  map  and  plan  measure. 

Stage -Area  and  Stage - Pis  charge  Felat ionshlps 

For  the  out f. low  section  of  each  of  the  sub'.’atersiieds ,  stage- area  and 


stage-discharge  curveo  must  be  supplied  to  the  PF3^.  To 
nation  it  mlis  necessary  to  obtain  channel  cross -sect ions 


tw ent y “G ix  local ions 


get  this  inf or- 
in  the  fie. If?  at 


iTom  the  cross-section  siurvey 
by  assuming  that  the  water  surface 
The  stage -discharge  ciurvs  for 


notes,  stage -area  curves  vaerc  obtained 
is  always  laterally  horizontal, 
eacli  section  was  ca3.cu.latcd  by  the 


kianning  equal  ion : 
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0  ar'’-667  cO.5 


n 


(33) 


v.’here  n  is  the  I-'snjiiDs  rou^hnass  coeffic:' ent ^  11  the  hydraiilic  radius ^  and 
S  the  friction  slope  at  that  section* 

Friction  SJ.o;u.2^i  Friction  si.opa  is  fallen  to  he  equal  to  the  channel 
dope  on  Fuck  Creek  for  several  reasons.  .  ... 

First  of  ail,  Henderson  (I966)  states  that  because  of  the  steep 
S-lopes  usually  present  in  the  upper  catchment  rep:ions  vhere  the  runoff 
prohlen  exists,  bed  slope  will  be  the  only  siynif leant  slope  teni  and 
the  discharge  viil  be.  a  function  of  depth  alone*  Bed  slope  can  he  con¬ 
sidered  the  significant  cloi)e  factor  for  slopes  greater  than  10  feet  per 


mi  Ic 


Furthermore j  the  cross-sections  on  Buck  Creek  arc  taken  at  an  average 


of  tuo  miles  apart*  Also,  there  are  many  variations  in  cross-section 
shape  hotveen  any  t\-;o  measured  cross-sect  ions .  Therefo^’e,  the  actual 
friction  slope  could  he  expected  to  deviate  from  the  average  bed  slope 
throughout  the  channel  reach.  The  difference  betueen  the  average  friction 
slope  and  the  average  bad  slope,  hovrever,  v.’ould  be  expected  to  be  vei’v 


Trans  it  ion ;  Tiis  stage-discharge  curves  used  are  a  corabina- 
tion  of  curves  based  on  different  assumptions  *  The  channel  slope  used  for 
lower  flows,  where  the  flow  is  confined  to  the  min  channel,  is  the  slope 
of  the  channel  bed*  P’or  very  large,  flows  covering  the  flood  plain  and 
f 2 owing  at  depths  greater  than  the  vegetation  height,  the  effective  slope 
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is  that  of  the  flood  pj.atn„  imxoortaiice  of  this  consideration  of 

slope  is  ax^parent  v;hen  the  ratio  of  meander  length  to  flood  X'lain  length 
Is  high;,  as  it  is  on  lov;er  portions  'of  Duck  Creek  Mhero  it  reaches  a 
value  of  tyo« 

A  transition  zone  exists  in  xhich  the  effective  slope  varies  ba“ 
tween  the  channel  bed  and  flood  plain  slox^as  (shown  in  Figure  9)®  This 
transition  would  be  difficult  to  simulate  without  extensive  knovLledge  of 
the  w’atershedo  Therefore;  in  this  study  the  stage  "'discharge  curve  is  a 
function  of  the  channe-l  bed  s.lox)e  until  approximately  one  foot  of  water 
is  flowing  on  the  flood  plain;  above  this  depth  the  flood  plain  slope  is 
cons idered  effect ive  c 

thnning ' s  Roughness  Coefficient ;  G'lie  values  of  Manning's  roughness 
coefficient  "n'‘  were  selected  using  a  number  of  photo^^rajfns  of  channel  , 
sections  for  which  "n”  had  been  calculated  by  \M  Th  Chov?  (.195Q)»  These 


v.’ere  compared  with  of  Duck  Creek  to  determine  the  n  value 

used  for  the  various  channel  reaches  of  Duck  Creek.. 

An  attempt  was  made  to  check  this  determinatioii  vi.th  a  set  of  U.  S, 
Geological  Survey  (U3G3)  stereoscopic  color  slides  that  have  been  used 
for  roughness  coefficient  detcimiinat ions  for  that  organization.  Unfor¬ 
tunately;  the  set  of  slides  viewed  did  not  include  exaniplas  similar  to 
those  found  on  Duck  Creek. 

One  further  comparison  was  made  in  an  attempt  to  evaluate  the 
effective  roughness  of  Duck  Creek.  This  method  was  developed  by  Cowan 
(3.95^0  can  be  siuruarized  with  the  following  rc.lat ionship : 
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Figure  9^  Typical  Plan  Viev  and  C;i 
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n  =  (hq  +  nj^  +  Iig  -I-  nj  1-  )nj  (34 ) 

vhere  n  is  the  desired  roughness  coefficient  to  be  used  in  ^hnning's 
formula,  uq,  n^,  n2,  n^,  and  n^^  are  the  fractional  values  of  n  deter¬ 
mined  respectively  by  type  of  channel  bed  and  side  slope  material,  effect 
of  surface  irregularities,  variation  in  shape  and  size  of  cross-section, 
effect  of  obstruction,  and  vegetation  conditions,  n^  is  a  multiplying 
correction  factor  and  is  a  function  of  the  ratio  of  meander  length  to 
flood  plain  length  for  any  given  reach  of  channel.  The  roughness  coef¬ 
ficient  given  by  this  method  is  larger  than  those  determined  by  photo¬ 
graph  comparisons. 

The  problem  connected  with  photograph  viewing  for  ”n"  determinations 
is  that  effect  of  meandering  cannot  generally  be  determined.  If  one  is 
using  the  channel  meander  length  on  slope  determinations  then  the  "n" 
value  found  is  probably  correct.  However,  if  the  flood  plain  length  is 
being  used  (which  is  the  case  for  overbank  flood  flows)  "n"  must  be  ad¬ 
justed  for  the  degree  of  meandering. 

For  the  in-channel  flow  condition,  "n"  was  determined  as  discussed 
above.  The  meander  channel  slope  was  used  for  this  condition.  However, 
for  the  overbank  flow  condition,  "n"  on  the  flood  plain  was  determined  by 
previously  explained  methods  and  was  higher  than  the  channel  "n"  used  for 
the  in-channel  condition.  The  slope  for  the  overbank  flow  is,  of  course, 
a  steeper  one.  "n"  was  taken  to  be  O.OO5  higher  for  the  main  channel 
portion  of  the  cross-section  than  for  the  overbank  flow  case.  This  is 
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done  to  account  for  the  fact  that  much  of  the  flou  In  the  area  of  the 
channel  is  actually  following  the  channels  Ihxs  ia  mox=e  convenient  than 
aesuminc^  that  the  channel  portion  of ■  the  flou  is  acting  on  a  milder  s.lopcc 

Eydranlic  P.acl  iu.s  i  Often  the  croGs-section  of  a  natural  strccam  is  so 
irreguj.ar  that  the  thnning  equation  cannot  be  conveniently  a]ppliec1  to  the 
croGB-section  as  a  whole.  For  the  water  surfecce  elevations  where  the  flow 
is  in  transition  betvjeen  channel  and  flooci  plain  flow,  •  neither  the-  wetted 
perimeter  P  nor  the  hydraulic  radius  R  accurately  represents  the  actual 
condition.  An  alternate  way  to  express  this  idea  is  to  say  that  the  ex-“ 
X)onent  on  R  does  not  equal  O.667  in  the  transition  zone.  This  is  gs-= 
pecially  true  if  the  main  channel  is  deep  and  the  flood  plain  is  qiii-te 
flato  To  alJ.evlato  this  problem,  the  cross-sectional  area  can  be  bro’uen 
up  into  several  portions  riebing  the  port. ion  divisions  at  the  sharp  breahs 
in  the  section. 

Coimputer  Solut ion  of  Manning ' s  Eguat ion ;  A  computer  program  was 
written  by  the  author  in  the  Fortran  II  3.anguage  for  the  IBM  l620  com¬ 
puter  (Appendix  B)  to  soirve  thnning's  Eqv-aticn  for  Q  at  a  given  cross- 
section  for  all  depths  (by  one -foot  increments)  between  spsecified  limits, 

Ifhere  the  cross-section  is  irregiular,  program  fircst  computes 
area  and  wetted  perimeter  for  each  portion  of-  the  cross-section.  The 
stage-disebargo  information  is  then  printed  out  for  that  portion.  The 
same  ca3.cu3.ations  arc  irade  and  printed  out  for  all  the  portions  of  a 
paj-ticular  cross-section.  Finally,  the  stage-discharge  values  for  all 
portions  ai'e  suirmcd  up  and  printed  out  as  the  stage-discharge  relationship 


valid  for  the  "rfliole  cro3G“f:action.  Tnis  GUFirary  information  is  then 
used  as  input  to  the  PFPo 

Tasting  the  Rellai.'bility  of  Manning’s  Equation:  The  reliahilit v”  of 


ihrining 'fi  equation  and  of  the  computer  prograifi  v?as  tested  by  comparin 


the  stage “discharge  rating  curves  uith  those  uhich  have  been  developed 
by  the  IJSG13  for  the  East  Pork  (Figure  10)  of  IXick  Creek  and  for  Duck 
Croak  (Figure  11)*  The  rating  curve  for  the  East  Fork  v’as  developed  for 
the  same  location  (previously  described)  used  in  this  study*  The  rating 
curve  for  Duck  Creek  as  developed  by  the  USG3  is  based  on  their  staff 
gage  located  at  the  bridge  v?here  Duck  Creek  flows  under  the  county  read 
between  sections  11  and  lb;  Tl'/N;  Rb6E.  The  Montana  State  Uaiversity 
stage  recorder  was  originalJ,y  placed  at  this  same  location  hut  was  moved 
ill  June;  1964  to  a  point  330  feet  north  (downotrean)  of  the  county  road 
centerline,  Tiie  cl.d  rating  curve  does  not  now  app3.y  to  the  new  location 
of  the  ISU  stage  recorder  and  since  this  paper  inc.lucles  a  study  of  the 
runoff  event  of  June  I6;  19^5  ^  i"®"’  rating  curve  had  to  be  developed* 

Flow  has  bean  iiisasurcd  on  three  occasiono  by  the  USGS  since  the  re¬ 
corder  on 'Duck  Creek  was  moved,  Ti^^o  of  these  measurements  were  made  in 
the  STiring  when  the  flow  was  over  an  iced  channel.  Therefore;  an  accurate 
rating  curve  cannot  be  drawn  based  on  actual  measurements  *  3'.t  was  do- 

terininad  that  the  Manning  equation  gave  a  sufficiently  accurate  rating 
curve.  This  was  done  by  comparing  the  actual  and  csjculated  rating  curves 
at  the  recorder  on  the  East  Fork*  At  this  location.;  the  rating  curve 
given  by  Manning's  equation  fell  near3.y  on  ton  of  the  best  visua3..ly  fitted 
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Figure  10:  East  Fork  of  Duck  Creek  Rating  Curve 
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line  through  polntB  maasured  hy  tlie  USG3.  Since  this  is  true?;,  it  my 
he  assumed  that  the  rating  cm've  for  Duck  Creek  calculated  by  the  i'anniDg 
relationship  should  be  adequate  for  pui’poses  of  this  study ^ 

As  alrea.dy  stated tuo  different  slopes  have  bean  utilized  for 
developing  rating  curves  using  ilanning’s  equation,  depending  on  \Jhether 
the  flcc’  is  inainl:/  in-channcl  or  on  the  flood  plain.  In  the  foriLer  case, 
the  length  of  the  niain  channel  centerllDG  or  the  meander  length  was  used 
to  deteri'iine  the  bed  slope.  In  the  latter  case,  houeve'^,  the  length  of 
the  flood  plain  v;as  used.  This  length  is  of  course  shorter  than  the  riaaudo 
length  giving  the  f.loodplain  a,  greater  slope  than  the  nain  channel.  In 
the  case  of  Duck  Creek,  the  average  ratio  of  meander  to  floodplain  lengoh 
is  about  1.89. 

Ci.u:ve  Nurabers  and  Antecedent  Mo  let  lire  Condition 

As  shoun  in  Chapter  III,  the  soil-cover  complex  number  (CN)  plays  a 
vital  role  in  the  SCS  method  of  determining  runoff.  A  composite  Cil  must 
be  calculated  for  a  watershed  if  more  than  one  soil  type  or  cover  con¬ 
dition  exists.  The  curve  numbers  used  by  the  SCS  for  various  combinations 
of  soil  and  cover  conditions  are  given  in  Appendix  C. 

A  composite  CN  vms  calculated  for  each  of  the  2T  subv/atersheds  on  the 
Duck  Creek  watershed.  The  method  used  to  compute  this  composite  is  shorn 
below  and  is  summarized  in  Table  3«  First  it  was  deteimrined  what  the  CN 
would  be  if  the  whole  subwaterslied  consisted  of  one  soil  t;)n:)s.  Tliis  was 
done  on  a  v;eighted  basis  according  to  the  percentages  of  acres  falling  in 
the  various  cover  classes.  This  was  done  for  each  soli.  tyx:e  (types  B  and 
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(3) 

(^) 

(5) 

(6) 

(7) 

Soil 

Cover- 

Curve  Kuaiber 

f.  of 

Column 

of  area 

Type 

Condition 

from  Appendix  C 

total  area 

(3)  >:  (!0 

of  soil  type 

(5)  X  (6) 

B 

Strip 
farmi  ng 

80.5^:- 

21c  4 

17.2 

lic'.uge 

69 

73. 6 

94.2 

Cl!  if  all  tn-^e 

B  soil  c . . . 

.  71.4 

29c3 

20..  9 

n 

Strip 

farming 

’87^^ 

21.4 

18.6 

Range 

79 

78.6 

62.0 

Cl’J  if  al],  type 

C  SOlJ.  ecct 

.  d0c6 

70.7 

97.0  1 

Coraposite  CN  (rounded  to  nearest  uholo 

number )  s 
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Calculation  of  Composite  Curve  iNumbar 
on  a  Sample  Subvatersbed 


*An  average  of  CN'c  for  fallou  and  small  grain  (straight 
hydrologic  condition)  cover  condition  classes. 
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C  v'ere  the  onjy  tj-yes  occurri;?g  on  Duel;  Creels). 

'j-he  weighted  CN'S  obtained  for  each  coil  type  individually  are 
then  again  v/eighted  with  regard  to  the  percentages  of  area  within  each 
subwaterched  that  fall  into  the  various  .soi.l  groups.  Tliis  method  of 
figuring  a  composite  CD  is  based  on  the  assumption  that  the  various  soil 
type  percentages  used  occur  in  those  same  percentages  on  each  of  various, 
cover  condition  groups  used.  If  this  is  not, true,  significant  errors  in 
CN  calculation  couJrl  occur.  It  is  felt  that  such  an  error  would  not  be 
important  in  the  case  of  Duck  Creek  since  in  nearly  all  subwatersheds 
one  soil  type  or  cover  condition  was  predominant.  Vlheii  this  is  not  true, 
it  j.s  necessary  to  calculate  the  v.’eighted  CN  for  each  subarea  of  eacl- 
subwatershed,  where u'pon  the  breakdov/n  of  soil  types  for  each  subarea  would 
have  to  be  determined.,  This  detaii  was  not  considered  necessary  on  Duck 
Creek . 

Tha  determination  of  CN  from  Appendix  C,  which  is  based  on  an  average 
antecedent  liioistvare  condition  (AMC),  is  input  to  tlia  PFP.  Also,  the 
actual  AKC  is  determined,  as  described  in  Charter  III,  and  read  into  tlic 
PFPe  The  adjustment  of  CK  is  internal  to  the  PFP«  For  the  storm  of  June 
l6,  1965  on  Duck  Creek,  the  AFC  v/as  determined  to  be  dry  (Type  I), 

Time  of  ConcerrLrat  ion 

T'l-e  time  of  concentrat ion  for  any  givtui  charnel  roach  is  taken  to 
be  the  time  of  travel  for  full  bank  flow  condltionG.  The  velocity  for 
full  hank  flow  is  determined  by  Idimi.ing’s  equation  at  both  ends  of  the 
channel  reach.  Travel  time  is  then  given  by  the  following  relationship: 
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Vj  +  V 


(35) 


I,  is  the  reach  length  in  feet;  Vj  and  Vp  are  the  -  velocities  at  the  in¬ 
flow  and  outflo'7  sections  of  the  channel  reach  and  t  is  the  travel  tiLie 
tbrongh  the  roach.  In  the  case  v/hero  the  roach  in  qiAsstion  is  between: 
a  st-ream  cross  section  and  the  boundary  of  the  v/atershed;  velocity  is 
ta.hen  to  be  '^ero  at  the  boundary.  These  times  are  then  added  down  through 
the  v?atershed  to  find  the  watershed  lag  time  (time  between  centers  of  liUiss 
of  the  rainstorm  and  the  outflow  hydrograph)  for  the  watershed ,  hTiCn  this 
calculation  does  not  match  the  actual  lag  tins  ar3  observed  from  outflow 
hydrogrepbS;  each  reach  travel  time  must  be  adjusted  by  proportional 
multiplication  so  that  the  actual  lag  time  is  obtained.  The  time  of  con¬ 
centration  for  each  individual  subwatershed  is  then  taken  to  be  a  poirtion 
of  the  stream  full-bank  fJ.ow  time  plus  overland  and  tributary  channel 
flow  times.  Tlie  combination  of  main  stream  end  side  channel  travel  times 
that  is  longest,  is  taken  to  be  the  time  of  concentration  for  the  sub¬ 
watershed  . 

The  tributary  channel  travel  times  v;ere  based  on  the  velocity  of 
two  or  three  feet  per  second  at  the  mouth  and  zero  at  the  head.  These 
two  velocities  and  the  length  of  the  tributary  chanrel  are  then  sufficient 
to  give  the  time  of  travel  as  given  by  equation  (3?)-  This  is  an  approxi¬ 
mate  method  but  sufficiently  accurate  on  the  Duck  Creek  .subwatersheds  be¬ 
cause  they  are  characteristically  loiig  and  narroAiU 


The  tributary  channel 
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•limsG  of  travGl  were  short  compared  to  the  main  channel  travel  tiraec. 

For  shorteih  v;lder  subv?atersbeds  the  deterrai nation  of  time  of  travel 
and  time  of  concentrat ion  v/ould  need  -more  scrutiny,  for  the  side  channel 
times  mould  constitute  larger  percentages  of  t-he  total  times . 

Storm  fata 

Both  actual  and  design  storms  can  be  handled  by  the  I'FPo  The  rain- 
fall  amount  is  always  i*ead  in  as  either  the  actual  amount  from  some  storm 
or  the  design  rairh'all  amount  as  dictated  by  c-oine  design  criteria*  The 
pattern  of  an  act\ial  storm  can  be  taken  into  accvrunt*  For  instance,  v?hen 
the  uatershed  is  broken  down  into  several  subwatershedo  a  different  rain¬ 
fall  amount  can  be  designated  for  each  v/atershed*  Ojempora.l  variations  can 
also  be  accounted  for.  in  actual  storms  *  Tne  starting  time  of  the  rain¬ 
fall  can  be  different  for  each  Guhwatershed* 

Balnfall  patterns  can  be  read  into  the  program  as  cumuJative  dimen- 
slonless  graphs  or  graphs  vritli  the  actua],  time  and  rainfall  amounts  as 
abscls.sa  and  ordinate,  respectively*  As  many  as  nine  different  patterns 
can  be  used  in  any  one  run  of  the  computer  prograra. 

The  amount  of  rainfall  falling  on  each  of  the  ?li-  subwaterslicds  on 
Duck  Creek  was  determined  by  the  Thiessen  Method,  The  Thiessen  polyg^cn 
was  drawn  for  the  seven  rain  gages  on  or  near  the  vvaterslied,  VJith  this 
method  the  rainfall  over  the  entire  Thiessen  svibai'sa  is  considered  to  be 
uniform  and  equal  in  value  to  the  araount  recorded  at  the  gage  wiiic’ii  gen- 
erally  falls  in  the  center  of  tlie  polygon  if  the  gages  are  equally  spaced. 
As  V70u.ld  be  expected,  some  of  the  suhwatersheds  fall  in  two  or  three  of 
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the  Thiessen.  polygons in  \which  case  the  raii-vTall  amount  is  a  v’cighteG 
amount  based  on  the  percentages  of  area  falling  in  each  of  the  po.lygons 
Fortunately,  the  edges  of  tbs  pcnlygcns  fall  close  to  the  i-idgo  be- 


twaan  the  East  and  VJefit  Folks  of  Duck  Creek:  this  eliminates  the  large 
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criticism  of  the  Thiessen  method  uhich  says  that  the  method  does  not  take 
into  account  geographic  features c 

Storm  movemant  was  taken  into  account  by  varying  the  starting  time 
for  each  of  the  subwatersheds «  The  speed  of  the  storm  was  determined  by 
exairdnation  of  the  records  of  the  three  continuous  recording  rain  gages. 
The  time  for  movement  of  the  storm  across  the  watershed  was  found  to  be 
about  3.5  hours  for  the  storm  of  June  I6,  1965e  The  startin.g  times  for 
the  various  subc’atershcds  were  then  determined  to  the  nearest  0.1  hour  by 
visually  determining  the  relative  positions  of  their  centroids. 

Summary 


The  amount  of  and  kinds  of  data  supplied  to  the  PFP  vary  frora  water-” 
shed  to  watenrshed.  Not  all  types  of  data  thc,t  could  be  processed  by  the 
program  need  be  used  in  every  vmtershed  model  developed  for  the  computer 
solutioru  For  example,  although  there  are  a  number  of  small  stock  pond 
reservoirs  on  Puck  Creek  watershed,  their  effect  is  felt  to  be  minor,  and 
therefore  the  model  does  not  include  any  reservoir  inuonnat ion.  The  data 
supplied  to  PFP  has  been  discussed  in  this  chapter.  With  this  information, 
the  PFP  is  capable  of  caJxu.lating  the  runoff  hydrograph  from  Duck  Creek 
based  on  the  usual  SCS  hydrologic  formulas  and  proce  11,1:1x6 s . 
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Chax)tor  VI 


ACIUAL  AND  liVTOTHETICA L  S^ODI'D 
The  actual  storm  of  June  l6,  I965  on  Duet  Creek  nas  used  to 
adjust  a  rainfall-- runoff  model  of  the  ^vaoershed»  'ifL&n  hyx^othetical 
storms  based  on  \L  B.  Weather  }3ujreau  frequency  data  v?era  routed  through 
the  iticdal  to  rn-ediet  peak  discharges  c 
Actual  Storm  of  June  196^ 

The  actual  storm  used  for  ir,odel  adjustment  on  Duck  Creek  has 
heretofore  been  cualled  the  storm  of  June  l6^  19^9 ®  dhs  reGultiug  peak 
dischairge  actually  occurred  on  this  dato^  although  the  stoxan  in  question 
started  at  about  6  a.m*  on  June  I3  and  lasted  unto’ 1  about  amn,,  June  16 
fox’  a  total  durattion  of  70  hours.  These  times  are  auex'age  times  as  the 
storm  had  slightly  vax-ying  dm’atioiis  as  it  ixoved  across  the  vjatershed. 

The  majority  of  rainfall  occurred  di’ring  three  main  birrsts  ap'proxi" 
mately  a  day  apart.  The  duration  of  the  bursts  uere  approximately  3;  1; 
and  6  houi'3, 

Figux’c  12  shous  the  area.!  distribution  of  the  total  amount  of  rain'* 


fall  in  the  foriri  of  an  Isolryetal  map.  This  is  based  on  the  recox'ds  from 
a  total  of  sevcii  x^aln  gages  on  or  near  the  v’atershed. 

The  storm  moued  generally  in  a  northwesterly  direction  down  the 
watershed’s  main  channel,,  An  indication  of  the  storm  movement  was  ob- 
tained  by  study ii'g  j'ecords  of  time  versus  cumulative  rainfall  amount  by 
poresntago  for  two  of  the  thx'ee  continuous  recording  x.-’a'In  gages.  These 
indicated  that  the  differonce  i)i  time  at  the  two  I’ccorder  points  at  which 


the  sam.3  percentages  of  total  sto?,'m  amount  had  c<cciirred  was  about  tv?o 
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hours  ^  This  Mas  apprcximtoly  truo  for  all  three  bursts  so  that  one 
storm  pattern  was  assumed  to  be  adequate  for  all  of  the  subwatersheds . 
Assuiaing  that  the  storm  traveled  at  a  constant  rate  of  speedy  the  feO 
hour  difference  in  tine  of  stona  beginnlrg  hetucon  the  continuous  re^ 
cording  points  bcccmes  3,5  houx-s  from  the  areal  centroids  of  the  two 
most  widely  sex^aratod  subwatersaedso 
Hypothatical  Stor-ms 

Storm  durations  of  12  ai-d  2k  hours  were  chosen  arbitrarJJ^y  for 
the  purpose  of  using  the  SC3  method  to  predict  peal:  discharges.,  rain¬ 
fall  amOi-vUits  routed  for  both  durations  were  j,n  0c5  inch  inci’cm.ents  and 
over  a  ra.nge  such  that  5^  -iO.?  25 >  50;  and  100-=  year  frequency  storms  woul.d 
be  included..  These  amounts  are  those  predicted  by  the  heather  B’uoeaai  in 
their  Technical  Paper  IIo*  1|0  as  X'^’epared  by  David  liershfiei^d  (1961).= 
Therefore;  for  12-luOiir  storms ;  precipitation  P  of  2.0;  2^5^-  3«0;  3«5^ 
and  4.0  inches;  and  24-hoiir  sto?cms  v/lth  a  P  of  2.0,  2.5^  3*0;  3e5; 
and  4c 5  inches  were  routed. 

The  time  of  storm  be,gimiln3  was  tahen  to  be  constani:  on  all  sub- 
v?at6:r.sheds  since  there  is  little  information  on  the  direction  of  pre¬ 
dominant  storm  movement  on  Duck  Creek,  If  a  consistent  pattern  of  storm 
movement  was  known  in  the  case  of  a  x^a^'l^icular  watorshed;  than  some  varied 
time  of  storm  beginning  could  bo  assumed;  even  for  hxqvothetical  storms. 


AlsO;  the  total  g-mount  of  rainfall  must  be  a.GSi’rad  to  be  corstair 
on  all  .subwatersheds  unless  there  is  infoiTiat ion  indicating  that  the 


storm  center  for  most  knou’n  actual  stoma  is  consistenbly  located  in  one 
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portion  of  the  naterohed*  Suoh  inf oiniiat ioii  is  not  kncnn  for  f/rck  Creek 
Therefore,  rainfai.!  auiQunt.s  vere  considered  to  be  spatially  invariazit., 
The  effect  of  the  location  of  storm  centers  could  be  studied,  hoi; ever, 
by  centering-  storms  in  various  parts  of  the  watershed  and  comparing  the 
calculated  outflci:  hydrograph  results ^ 

Actual  Runoff  Hydrugi’aph 


The  actual  runoff  liydrograph  for  the  c-toriii  of  June  16,  1905  is 
shoim  in  Figure  Id®  This  hydrograph  had  been  developed  by  previous  in¬ 
vestigators  based  on  the  best  determination  of  tlue  rating  cu.rve  in  effecl 
at  the  new  streani  stage  recorder  location®  The  volume  of  ruinoff  associ¬ 
ated  with  this  curve  is  0.^3  inches  for  the  9^  sejuare  mile  Buck  Creek 
watershed . 

After  this  study  i-as  \/ell  under  way  it  v?as  discovered  that  two 
errors  had  been  coifim.ltted  in  the  reduction  of  the  continuous  stage  re¬ 
corder  data,  Tlia  first,  mafle  by  Etudv-riuS  working  on  the  Drainage  Curre- 

J.ation  Research  Project,  was  that  the  volume  previously  calculated  was 
in  error.  The  I'esu-lting  yield  should  ho^ve  been  0,28  instead  of  the  0,b3 
inches,  liio  eacoi'd  error  was  made  by  the  U.  S,  Geological  Surrey  In 
thsir  publications.  The  peak  stage  printed  for  the  June  I6,  I965  storm 
was  A, 7’'-  feet  v?hlch  gives  a  jjoab  of  8'{0  efs  based  on  their  rating  curve 

for  Buck  Creek o  The  vStage  value  should  have  been  2,83  feet  which  gives  a 

peak  of  3^+0  cfs  instead.  Since  the  original  ro.ting  curve  fur  the  now 


location  was  based  prirrarily  on  this  one  storm,  it  is  now  invalid, 
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ThG  author  has  clave  .loped'”'  another  .ratiug  curve  based  on  the  Lh’rAr.in^^ 
Equa'tiO);h  this  bain-^  just'J.f'ied  for  reasons  dcBcribed  in  Chapter  Y,  Onl 
tuo  points  are  kncuu  on  the  rating  curve.  One  is  the  stoi-’n  being  tc;::.ted 
and  the  other  is  the  elevation  or  stags  at  v?blch  the  flow  is  j'.ero,  A 
ra'bing  curve  uas  drawn  through  thecc  tuo  points  based  on  the  hrinnirig 
Equatlnr-c.  The  runoff  hydrograph  uas  again  dravn  based  on  the  new  rating 
curve ;  and  is  also  chovn  in  I''igiire  13.  Ths  ordina’oGs  of  this  plot  o’ers 
read  into  a  coup'ater  program  writc-en  by  Gary  Iruis  to  determine  the  yie.ldo 
The  yicl.d  non  b-3.1  ieveu  correct  for  the  June  l6^  1965  aterm  on  Duck  Creek 
is  0,35  inches, 

Coiaputer  Rjyn;£  Comp le ted 

A  number  of  rims  or  storm  routings  were  made  on  tvjo  dif'ferent  occa¬ 
sions,  TnsSG  runs  are  sur^imariced  in  Tab.le  4, 

Runs  1"*A  through  l-M  and  2 -A  tlu’ough  2-M  are  identical  except  for 
the  curve  numbers  (CH)  used.  The  nuLbor  .1  runs  rsake  use  of  a  set  of 
CR's  criginal.ly  tiiought  uO  be  those  in  effect  on  Duck  Cx’gs'k,  For  the 
number  2  runs,  each  CN  of  the  sot  uas  reduced  by  the  value  of  th.rce* 

The  A  and  B  runs  (both  numbers  1  and  2)  caused  the  June  l6,  19^5 
storni  to  be  routed  tVircugh  the  Diiclr  Crock  noclcl.  The  only  difference 
bstyo’cu  A  and  B  is  the  asc-umption  made  about  the  t'ima  of  the  storm’s 
bcginnlrgo 

Runs  C  through  M  arc  for  hypothetical  stormis.  Runs  C  to  H  arc 


•»”V/ith  assistance  freca  Gary  Irr.rls,  graduate  student  at  lioTilana  State 
University  working  tcuards  an  M.  S,  Dsyx'ee  in  Civil  Engineering. 
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Run 

of 

run 

T7^pe  of 
Storm 

i 

Amour y 
(in,) 

Barat ion 
(hours  ) 

Time  of 

St  orm ' s 
Beginning 

. 

Antecedent 

MoiGtiire 

Condition 

Curve 

ITuttBc  3’ 

CN 

1-A 

j.o/ 

Actual 

Va  r  i  c  1. 

70 

Varied 

E-y 

Griglrally 

E 

25/ 

Actual 

Varied 

70 

Constant 

Ticy 

doteimn.ncd 

C 

67 

24 

Constant 

Average 

from  AppsnCiix 

3) 

tlietl- 

*  0 

C  for  each 

E 

cal 

oo 

subvaterahed 

F 

3.0 

G 

H 

2,0 

I 

4,0 

12 

J 

3o 

■ 

1 

1 

K 

3*0 

i 

L 

2,5 

1 

Id 

2  VO 

< 

i 

) 

: 

2  -A 

11^ 

Actual 

Varied 

70 

Varied 

Dry 

CN  non  three  | 

B 

21/ 

Actual 

Varied 

70 

Constant 

3):i’y 

less  tlian  | 

C 

67 

Hypo- 

4 , 5 

24 

Constant 

Average 

those  for  runs 

1) 

theti- 

4,0 

l-7\  to  1-M 

E 

eal 

3.5 

F 

3*0 

G 

2.5 

11 

2,0 

I 

4,0 

12 

J 

3c5 

K 

3,0 

T. 

2,^ 

h| 

2,0 

1 

i 

Table  ^i- :  Sumnary  of  Corq')Ut-er  Euno  Completed, 
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24*»hoiu’  duraticn  ctorBS  vith  varl^abls  rainfall  amounts^  vhila 
are  12 hour  diu’ation  storms  with  variable  rainfall  amounts.  ■ 
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ESSULTS 


Tlie  data  available  on  the  storu  of  June  l6;  196^  on  the  Jlick 
Creek  watershed  v’as  used  to  set  up  the  coiaputer  uodel  for  analysis  by 
the  Soil  Conservcition  Service  •»  Project  Formulation,  Program  (SCS-PFP), 

I.Vo  sets  of  computer  runs  uere  mde  by  the  SCS  in  Ilya  btsville  ^  Maryland 
on  an  IBM  36o/40  electronic  coinputero 

file  initial  run  gave  a  total  watershed  yield  of  0*37  inches c  Ihi£ 
uas  higher  than  the  actual  yield  measured  at  the  stream  gage  v?hich  was 
0*35  inches o  Hovjevery  at  the  time  the  com.puter  runs  v/ere  lus^de  the  actual 
yield  was  tiiought  to  be  0,^!-3  inches  (due  to  xrreviouyly  explained  error  )*, 


liie  model  was  initially  set  up  by  evaluating  the  soil-cover  coi-> 
I)lex  numbers  (Cil)  thought  to  be  those  in  effect  on  Duck  Creek*  liiese 
are  based  on  the  chart  used  by  the  SCS  (see  Appendix  C  of  this  paper)® 
The  procedure  used  by  the  SCS  is  to  check  the  yield  for  the  originally 
determined  set  of  CN's.  Tlis  CiT's  are  then  adjusted  uuuil  the  calculated 
yield  equals  the  actual  yield*  This  is  considered,  to  bo  the  calculated 
runoff  hyclrograph  which  can  now  be  compared  to  an  actual  runoff  hydro- 
graph  if  it  is  available c 

In  an  effort  to  correct  the  yield  to  0*03  inches  (thought  at  the 
time  to  have  been  the  actua.1  yield)  the  author  reduced  the  CD's  for 
all  subwaters hods  by  a  value  of  three  and  the  second  set  of  computer 
runs  was  lindc*  Ti.ie  resulting  output  v.’as  hydrographs  having  yields  ex- 


a  e t  l.y  e  q  ua  1  t  o  0.43  i  nch  0  s  * 


7 


Ccvjip'dtcr  Resd-lts  j  Figiiro  3,4  Oiiir.rG 


actu.a.l  runoff  Rydrograph 


for  the  JmVe  3.6,  3.965  sto.rm;  points  -indi eating  the  comijutcid  hydregraph 
peaks ^  and.  tuo  hydreycaphs  vfirich  have  been  p03tu.l,ated  oa  the  basic  of 
the  best  evidence  availab.le^  The  actual  bydrograph  as  measured  has 
a  louer,  longer  profile  than  tlie  calcul.ated  hydrograph.s «.  The  cal- 
ciliated  b.ydrog;raphs  have  bipli  peaks  and  steeper  recession  limbs* 

The  four  Isolated  pointc  shown  on  Figure  34  are  the  highest  peak 
dischargcivS  calculated  by  the  computer  runs  l-fy  1~B.  2-A_j  and  2‘B*  ‘The 
first  runs^  each  with  a  yle3-d  cf  0*57  inches are  indicated  by  the  two 
highest  points  (l-A  a-nd  l-B)^  The  higrier  cf  these  (l-A)  is  for  the,  case 
where  the  time  of  storm  beginning  'i?ac  assumed  to  vary  from  one  cub- 
watershed  to  the  next*  Tiiis  liydrogra-ph  has  a  peak  discharge  of  I4l6 
cfs.  The  lower  of  the  two  (l-B)  with  a  discharge  of  3245  cfs  is  based 
on  the  assumption  that  the  starting  time  v?as  the  same  for  all  sifowater- 
sheds*  The  later  computer  runs  with  reduced  curve  number  va3,uc3  re¬ 


duced  the  yield  to  0*^13  inches,  and  the  iieaks  (2-A  and  2-13)  to  13.44  and 
1006  cfs  respect ive.ly*  Again  the  higher  peak  v;as  for  the  case  of  the 
variable  storm  beginning* 

The  two  hydrographs  siiown  in  their  coinp3.ete  forra  on  P'igure  l4  ho.ve 
3.'eG?i  postulated  by  the  author,  and  are  based  on.  the  rosu3,.ts  of  the  3u.st 
cornputer  runs*  Each  of  these  hydrograx:)hs  has  a  yield  of  0*35  inches, 
matching  the  yie3.d  of  the  actual  hyaregraph.  Their  shapes  are  identical 
to  those  of  the  3rdter  computer  runs.  Iff.’  lurblier  computer  runs  should  bf' 
made  with  lower  CN  values,  it  is  expected  that  the  chapes  of  the  hydi’O- 
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graij'iw  might  cliaDgo  slightly.  Comparison  of  the  complete  hyclrogrcaplis 
of  the  four  computer  oruno  vhich  v?ere  nm.ds  indicates  that  the  change  in 
shape  would  not  be  great,  however. 

As  can  be  seen,  there  is  a  large  discrepancy  between  tiic  calculated 
and  actual  hydregraphs  with  regard  to  peak.  The  SCS)  states  that  the 
methed  was  not  designed  to  raatch  the  hydrograph  in  shape  but  only  in 
peak  value.  However,  in  this  case,  the  calculated  peak  is  either  2,4l 
or  2/l'^'r  times  as  great  as  the  actual  peak  depending  on  which  assumption 
is  used  for  the  time  of  storm  beginning  on  each  subw’atc-rshedo  Ihe  rami" 
fications  of  these  results  will  be  discussed  in  Chapter  VIII, 
H?,rpothctica3.  Btcnmi  Evaluation 

One  purpose  of  this  study  v;as  to  predict  flood  peaks-  for  various 
frequencies,  Ihls  was  accomplished  by  routing  a  number  of  hypothetical 
storms  through  the  v;atershed  model  to  determine  the  peak  discliargvC  pre¬ 
dicted  by  the  SCS  methodc  The  final  results  of  this  portion  of  the  work 
is  suinmarized  in  Fig;ure  15,  Of  course  no  actual  storms  have  occurred  by 
which  these  discharges  can  be  compared, 

Figui’c  15  is  tentative  a.nd  should  not  bo  used  for  design  purposes , 

It  is  show’n  here  only  as  a  demonstration  of  hov’  the  GCS-PFP  results  of 
routing  of  hypothetical  storms  can  bo  used  In  conjunction  v<itb.  the  U,  S, 
Vfeather  Bureau  frequency  data  to  x^redict  peak  discharges  from  small  watex 
EliedSc  The  two  curves  mai’ked  WB  (for  V/eather  Bureau)  are  plots  of  ex- 
XDseted  precipitation  P  versus  the  return  xpariod  In  years  (abscissa  scale 
at  the  bottom  of  the  graph)  for  the  storm  of  magnitude  P  for  both  24- 
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and  12-hour  storip.s  as  indicatoc.e  Tha  other  two  curves  liiarked  3>C  (for 
Duck  Creak)  suoi:  the  plots  of  P  versus  peak  dlscharije  qp  (abscissa 
scale  at  the  top  of  the  graph)  for  the  Duck  Creek  vjatershed  as  datei-- 
mined  by  the  computer  model  using  SCS  hydrologic  design  criteria  for 
the  same  tuo  stern  d’Vi’ationSc*  These  two  curves  rex)recent  the  modal  as 
set  up  for  the  later  computer  runs  (not  as  yet  adjusted  as  earlier  ex¬ 
plained  and  thus  not  suitahie  for  design  use)*  If  additional  adjust?oent 
of  the  juodal  could  be  mado^  the  same  procedure  as  described  here  uould 
be  appl5.cabla  in  setting  up  a  similar  set  of  curves  for  use  in  design.. 

If  Figure  15  v/cro  to  be  used  for  design^  the  follov;xng  procedure 
should  be  fol.loued  to  properly  predict  the  peak  dischai-ge  for  a  sto.rm. 
Generally,  strucuxu'es  on  a  v.’atershed  arc  hydraulically  designed  to  func¬ 
tion  properly  for  some  peak  discharge  caused  by  a  storm  of  a  particu.lar 
frequency  or  retxu’n  period*  Suppose  for  an  example,  as  shoun  on  Figure 
15^  ve  plan  to  design  a  culvert  that  v?ill  handle  a  5^^'"yGar  frequency 
storm  of  a  dux-ation  of  24  hours.  The  procedure  is  to  start  at  the  bottom 
of  Figure  I5  v?lth  the  50“year  storm  and  proceed  vex'tically  along  the  line 
shoun  until  the  WE-24  lioui'  curve  is  net.  Then  if  the  ordinate  P  is  I’ead, 
Me  find  that  the  Weather  Bureau  predicts  that  the  precipitation  is  3‘T7 
inches  for  a  yO-year  -  24-’hour  duration  storm.  If  hov.’ever,  v?c  follox  a 
horizontal  line  (to  the  right  in  this  examp.le  )  until  uo  meet  the  dis¬ 
charge  for  a  24-horu‘  storm  foxy  Duck  Ci'cek,  and  then  proceed  vertically, 
v/e  find  that  Op  r  3200  cfs. 

On.  the  other  hand,  -jf  the  peak  discharge  for  a  runoff  event  had  been 


measured  on  Duck  Creek 
cribed  in  the  previous 


;  the  exact  opi)Osite  procedure  to  that  des- 
paragraph^  could  be  followed  to  determine  the 


frequency  or  return  period  of  that 
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Chapter  VIII 


Discussion 

The  investigation  described  in  the  preceding  chapters  has  been 
of  two  types.  First,  a  model  of  Duck  Creek  watershed  was  developed 
and  then  tested  by  using,  as  input  data  to  the  computer  program,  in¬ 
formation  describing  an  actual  storm  which  occurred  on  the  watershed 
on  June  l6,  196^,  The  resulting  synthetic  hydrograph  generated  by  the 
program  was  compared  with  the  actual  hydrograph  recorded  at  the  mouth  of 
the  watershed.  The  model  was  adjusted,  after  initial  comparison,  and  the 
program  re-run.  Second,  the  Duck  Creek  model  was  used  to  predict  the 
hydrograph  from  certain  design-type  storms  based  on  V/eather  Bureau  storm 
frequency  data. 

The  Storm  of  June  l6,  196^ 

The  results  from  the  computer  runs  for  the  June  I6,  I965  storm  do 
not  compare  very  closely  with  the  actual  hydrograph  recorded  for  the 
event.  The  discrepancy  is  thought  to  be  mainly  the  result  of  imper¬ 
fections  in  the  model,  but  there  is,  however,  a  possibility  that  the 
actual  hydrograph  is  not  entirely  correct.  The  peak  discharge  that 
actually  occurred  is  felt  to  be  accurately  known,  because  it  was  de¬ 
termined  from  the  USGS  rating  curve  for  their  crest -stage  gage.  This 
rating  curve  is  based  on  ten  current  mater  measurements  over  a  period 
of  at  least  ten  years.  Except  for  the  peak  discharge,  the  rest  of  the 
actual  hydrograph  was  determined  by  relating  depths  recorded  at  the 
water  level  recorder  to  discharges.  As  has  been  mentioned  earlier,  the 
water  level  recorder  is  located  several  hundred  feet  downstream  on  Duck 
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Creek  from  the  crest-stage  gage,  and  the  stage -discharge  relationships 
for  the  two  stations  are  quite  different.  The  best  rating  curve  avail¬ 
able  at  the  water  level  recorder  uses  running's  equation  and  is  drawn 
through  only  two  known  points. 

Possible  imperfections  in  the  model  seem  to  fall  in  one  of  the 
following  categories:  1)  preparation  of  the  input  data  and  2)  the  SCS 
method. 

Preparation  of  the  Input  Data :  Wnen  one  is  dealing  with  hy¬ 
drologic  problems  the  situation  never  occurs  that  there  is  an  over 
sufficiency  of  data  available.  There  is  always  a  limitation  on  the 
time -dependent  data  available,  and  an  economic  limitation  on  the  thor¬ 
oughness  with  which  the  watershed  parameters  can  be  obtained.  A  rather 
significant  amount  of  watershed  data  were  obtained  for  the  study  re¬ 
ported  herein,  including  the  results  of  a  medium  intensity  soil  survey, 
numerous  infiltration  tests,  a  complete  stadia-level  survey  of  the 
principal  streams,  and  more  than  20  channel  cross-section  surveys.  It 
would  seem,  that  the  watershed  parameters  were  delineated,  therefore, 
more  completely  than  could  be  expected  in  the  typical  small  watershed 
design  investigation.  Results  of  the  study  show  that  even  this  pre¬ 
cision  may  not  be  adequate  to  insure  great  accuracy  in  duplication  of 
an  actual  hydrograph. 

The  way  in  which  various  watershed  characteristics  were  utilized 


will  be  discussed  below. 
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Tiie  cur-v'e  niiEiibers  (Ci'O  vjhich  are  in  the  SC3  program  to  charac-- 

teriza  coil  type  a/xl  land  uoe  treatnect  developed  from  J.arge  quan¬ 

tities  of  data  which  were  oui/ained  from  all  parts  of  tlie  United  States. 
Since  they  represent  s.verage.  conditions ^  and  ware  not  developed  speci- 
ficall5"  for  Montana^,  It  seems  likely  that  use  of  the  CU  table  in  Appendix 


C  may  be  responsible  for  some  cf  the  discrepaiicy  in  the  recnl.ti 


In 


particular^  it  seems  that  the  adjustment  which  is  madQ  for  dry  antecedent 
moistuxc  conditions  (type  I)  may  not  be  large  e.aough..  It  may  be  that 
for  Montana  w-atershedS;  the  adjusticent  bet-oeen  average  and  dry  conditicns 
should  be  larger  than  that  used  nationwide* 

The  CN's  used  in  the  first  set  of  computer  ru.ns  v/ere  adjusted  for 
the  second  set  of  runs  by  reducing  each  nirmber  by  a  value  of  three.  This 
reduction  theoretically  should  have  been  proportional  to  the  original  value 
of  each  Cil  for  each  subvatershed .  It  is  uncertain  what  the  effect  of  re“ 
ducing  in  equal  instead  of  proportional  amounts  might  bso 

There  are  indications  that  tiie  timing  of  the  runoff  in  the  model 
was  in  error.  This  seeiiis  to  be  especially  true  in  the  case  of  the  high- 
est  p>eak  shewn  on  Figure  id.  As  can  be  seen  from  the  main  portion  of  the 
actual  hydrograpli;,  lesser  peaks  occurred  both  before  and  after  the  rraxi- 
num  discharge.  The  calculated  hydrographs  seem  to  have  grouped  all  three 
peaks  5nto  one  and  consequent .ly  gave  too  hig)i  an  estimate  for  the  inaxi- 
mura  discharge.  Tiae  factor  that  controls  the  tir-rlng  on  the  model  is  the 
tirca  of  concentration  which  is  determined  for  each  subwatershed,  As 


TTientioned  earlier,  only  one  time  of  concentration  vcan  be  supplied  for 
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each  subv.’a'tershed^  and  thsrafcrc  no  variation  can  ba  assunicd  in  Ihic 
tirrn  for  differerri}  diGchar^ec,  The  time  of  concciitraticn  can  be  shewn 
to  vary  uith  rate  of  discharge  so  that  this  becoines  a  limiting  factor 
in  the  construction  of  the  modelc  Tne  times  of  concGiitrat ion  for  this 
study  v?cre  basc-d  on  velocities  as  calculated  by  the  fanning  relationship 
for  full  banl:  flcv-.  T’ne  acbual  event  in  19^5^  produced  little  or  no 
flooding;  it  seems  probable  therefore  that  the  channels  were  less  than 
bankfull  most  of  the  time.  The  true  velocities  uere  probably  less  than 
those  calculated,  and  the  times  cf  concentrat icn  J.arger  thran  those  cal¬ 
culated,  Ddiis  seems  -likely  to  account  for  the  fact  that  the  calculated 
peaks  occur  at  times  earlier  than  the  observed  peaks. 

8C3  Mat  hod :  The  SC3  method  v;as  developed  before  electronic  com¬ 
puters  could  be  utilised  to  solve  such  problams  and  therefore  had  to  be 
simply  constructed  so  that  it  could  be  solved  by  hand.  And  as  is  true  of 
all  methods  developed  for  predicting  peak  discharges,  the  SCS  method  is 
based  on  Lnny  assumptions.  Assumptions  are  made  either  because  5 1-  is 
thought  that  they  uill  not  greatly  affect  the  solution  or  because  they 
must  be  made  to  make  the  solution  economically  possible.  I-iiny  times  in 
solving  hydrolcgic  prcblems,  assumptions  must  be  mde  because  there  is  a 
scarcity  of  information  available even  though  it  is  recognized  that  these 
assumptions  rmy  significantly  affect  the  solution. 

The-  3C3  methed  can  be  thought  of  as  a  mathematical-graphic  i.iolel 
which  consists  of  a  number  of  formulae  and  graphs.  Although  some  flexi- 


bi.lity  has  been  bulJt  irlo  the  operation  of  the  PFP,  it  appears  that  the 
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£CS  nethcd  hss  not-  been  basically  altered  in  the  conversion  inlo  a  com¬ 
puter  nodel.  The  question  that  then  arises  is  whether  rhe  basic  r.ethod 
itself  cOvlIc.  now  be  tiade  more  flexible  for  computer  solutions. 

Points  in  the  model  that  seem  to  be  toe  rigid  aic:  1)  initial 
abstraction.  2)  leg  time,  3)  oynthotic  milt  hydrograph  shape,  4)  storm 
pattern,  and  p)  rcsemvoir  and  flood  routing  methods. 

The  initial  abstraction  is  assumed  to  be  equal  to  0.2  S  and  th? 
lag  tiir.e  to  be  equal  to  0.6  of  the  time  of  concentration.  The  coeffi¬ 
cients  of  0.2  ana  0.6  are  certainly  for  the  average  conditions  for  a 
number  of  v/atersheds  that  vere  studied  by  the  £C3.  The  rau  data  from 
which  these  coefficients-  were  determined  are  not  given  in  the  SCo's 
National  Engineering  Handbook  (N£H).  If  these  data  had  been  supplied,  it 
night  be  possible  to  choose  more  nearly  correct  factors  for  a  specific 
area. 


Also,  the  NEH  vdefines  only  one  synthetic  hydrograph  shape,  which  is 
intended  to  be  used  in  all  cases.  The  PF?,  heuever,  is  capable  of  absorb 
ing  any  unit  hydrograph  that  the  user  wishes  to  supply.  From  this,  it 
would  seen  that  the  SC3  method  as  outlined  in  the  NSH  is  more  restrictive 
than  is  necessary  nor  that  the  more  flexible  computer  programs  are  avail¬ 
able. 

Only  two  possible  stoimi  patterns  are  given  by  the  SC3  for  use  in 
testing  hj^pothetical  storms.  The  pattern  most  com-^only  used  (s/iovn  in 
Figure  5)  is  one  in  which  the  main  burst  of  the  storm,  comas  at  about  the 

.  The  other  has  the  rain  storm  biu’st 


midpoint  of  the  stom  in  time 


come 


earlier  in  storm.  Again,  the  PAP,  being  able  to  accoi;.modate  any 
storm  pattern  shape,  is  more  flexible  than  the  method  described  by  the 
SCS, 


The  HEH  suggests  several  nothcds  of  flood  routing  uhich  mriy  be  used. 
Only  two  of  these however,  have  been  incorporated  into  the  computer  pro¬ 
gram.  The  convex  method  is  the  only  procedure  for  channel  flood  routing 
available,  and  the  storage  “'indication  ir.atbod  is  the  on.ly  procedure  3'’or 
reservoir  flood  rouuing  available,  Tiiese  are  methods  wlricVi  are  easily 
adaptab3.e  to  manual  solution  and  perhaps  were  chosen  for  this  reason. 
Possibly,  other  methods  would  be  mor-e  realistic  and  as  easily  handled  by 
the  computer,  Ey  simply  adding  other  routing  subrcutlnes,  the  overall 
computer  model  could  be  made  more  flexible. 

Pcs Ign- Type  Stormc- 

An  example  of  design  peak  discharge  determination  based  on  frequency 
of  storm  was  given  in  Chapter  VII  and  needs  no  further  amplification  here 
It  might  be  pointed  cut,  however,  that  the  accuracy  with  which  the  peak 
discharge  of  a  hypothetical  stoT‘m  can  be  predicted  can,  be  no  better  than 
the  accuracy  of  the  determination  of  the  amount  of  rainfall  to  be  ex¬ 
pected  in  any  given  location  for  the  storm  frequency  chosen. 

The  hydrogi’apho  which  were  determined  for  design  stcrnvs  should  be 
considered  only  for  their  qualitative  mieaning  and  not  their  quantitative 
results.  The  reason  for  this,  of  course,  is  that  the  mwjdel  of  the  Wciter- 
shed  was  not  thorou’ghly  adjusted  in  the  sense  described  in  Chapter  III, 
Furthei’  adjustments  in  the  model  or  watershed  parameters  could  now  be 


.  ^  i  ■  .  . ;  '  '  i  Li 

■  1  ■  ^  y  •  Ti  ,  . 

.i  K  Hi 

V  f</j  Lp>A 
* 


1 


f 


■t 


•'i 


V 


V 


» 


A 


4 


4 

r 


i 

I 

\ 

I 


i 


*  * 

/ .  •  i 


\ 


■f/  ■■  '  .  I  j  ;•  I  »4'j 

.  fi;  ;.»* j  xij  Ji;  ^  . 

,  ^ 


’•  .  ‘  •  j 

■  >  '1  4  .  i"  ,  T 

L'  '*  ‘.io  '//  ’I'-TiJ  ;oi.  •-•  .•  6r,J*f 


. .'  -  ji' 


»  84.  ^ 


iiada  quite  readily ,  V/j/tli  a  satisfactorily  adjusted  inodel  available^  the 
hypothetical  storms  could  be  i'Ou.ted  through, th.e  uoclel  once  more,  and  de-- 
sign  hydrographs  «ou?.d  be  thereby  obtained.  As  further  lunoff  events 
occur  on  Duck  Creek.,  it  rray  be  conic  possible  to  determine  ho:.?  the  ruuoff 
pattern  changes  with  size  of  storm;  this  would  permit  modifications  to 
be  made  to  the  uatershaa  model,  which  could  then  again  be  subjected  to 
the  hypothetical  stoims,  and  still  more  definitive  resulbs  obtairud. 


itrogi'am  Use  by  Montana  State  Highway  bepartmenl'i 


‘fhe  resulbs  of  the  study  indicate  that  before  the  Iiontana  State 
Highway  Dspartinent  could  re3,iably  use  the  SCS  method  of  predicting  pro¬ 
bable  peak  discharge'  from  i-ontana  watersheds,  a,  long  range  study  and 
coLlection  of  data  would  first  be  rivccessary.  Before  the  BC3  method  could 
be  apxilied,  maps  locating  soil  typies  and  land  use  on  tbs  watershed,  or 
a  slngJ.e  ma];)  delineating  the  soil-cover  complex  niimbars  (Cll)  on  the  water¬ 
shed,  would  have  to  be  iHide  available,  land  use  rcajis  cou.ld  be  drawn 
largeJ.y  with  the  use  of  aerial  photographs  and  a  minimal  firsthand  know¬ 
ledge  of  the  locati.on  invo.lvcd,  HydroJ.ogic  soil  types  can  -present. ly  be 
determined  in  riiany  areas  in  Montana  from  existing  data  from  SCS  soil 
sm.‘veys,  Ilow’ever,  .large  sections  of  Monbana  have  not  yet  been  surveyed. 
Therefore,  a  comptlcte  sol-1  type  max.)  for  the  state  could  not  nov?  be  com¬ 
piled. 


For  predict  ion  of  a 


aesign 


dif 


charge 


from  a  particular  w'atersh 


additional  data  would  have  to  be  collected, 
channel  reach  .Tengthc  cou.ld  he  Hj-easured  from 


Areas  of  subwatersiieds  and 
aerial  xAiotographs  if  they 
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were  avai.lable  for  tlio  watershed  in  quest] oiio  Field  surveys  would  be 
necessary  to  deterraine  channel  cross  “-sections  *  Channel  bad  g1ox;g,s  would 
also  have  to  be  determined  in  the  fi-eld  unless  suitable  topcyiaphic  maps 
were  available  o  The  I-hnniDp;  roughness  coefficient  could  be  determined  by 
inspecting  the  area  or  by  viewing  photographs  of  the  channelo 

l-ontana  State  Universit^^  is  pres-ently  installing  a  third^gc aeration 
computer^  the  BD3  Sigma  J ^  which  will  be  caj>able  of  runxiinp;  the  current 
(1907)  version  of  tl}e  SCS'-FF?  with  a  few  rncdif ications  „  Unfortunately^ 
the  Sigma  7  v?aB  not  available  during  the  time  of  this  study ^  and  it  was 
not  feasible  to  make  -iuore  than  the  t'>;o  sets  of  runs  previously  described^ 
To  effectively  make  use  of  such  a  program^  the  hydrologist  should  have 
easy  access  to  the  cor/iouter  because  of  the  several  trials  v;hich  seem  to 


be  necessary  in  adjustment  and  use  of  a  modal* 

In  summary,  the  SCS  method  i*-equires  collection  of  a  roasonabl}'  large 
amount  of  data.  It  has  not  been  p)03Gible,  with  the  data  v?hic1i  v/ore  avail- 
able  for  this  study,  to  evaluate  the  SCS  method  to  determine  whether  such 
amounts  of  data  collection  are  justified* 
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Chapter  IX 


conclusions 


The  invsGt D.^at ion  reported  heroin  has  been  a  study  of  tho  SCS 
msthorl  of  predicting  peak  discharge  flow  lates  from  rain-caused  otorrr-s 
on  small  v/atersheds.  Data  for  one  actual  stc:rm  which  occurred  on  Duck 
Creek  in  Prairie  and  DcCone  Counties,  Montana,  together  with  several 
hypothetical  storms  for  the  sans  watershed;  were  processed  by  a  com¬ 
puterized  Eoluticn  of  the  GC3  rmthod.  The  resulting  sy.thetic  hydro- 
gj'aph  for  the  actual  storsi  was  compared  with  the  actual  hydrograph  whicli 
was  recorded  for  the  event* 

Findings  from  the  study  lead  to  the  follow’ ing  conclusions; 

1  »  The  SCS  nethod  is  a  systematic;  logically  organized  procedure 
v’hich  generates  a  synthetic  hydrograph  for  a  snail  watershed.  .It 
utilizes  data  characterizing  the  causative  rainstorm,  watershed  para¬ 
meters;  a.nd  accepted  flood  routing  techniques.  Ihe  roetbod  has  been 
effectively  progranjned  by  the  SCS  for  computer  solution. 


2  -  Comparison  of  tho  synthetic  hydrograph  generated  by  the  computer 
prograra  for  the  stern  on  June  l6;  I965  on  Duck  Creek  with  the  actual 


hydrograx'h  which  was  recorded  for  the  sane  event;  shows  considerable  dis¬ 
crepancy,  This  is  thought  to  be  X-^2.*imarily  caused  by  using  improper  value; 
for  \?atershed  paiameters  and  flood  routing  coefficients  as  input  to  the 
computer  program.  Alternative  factors  \;hich  rray  explain  part  of  tho  dis- 
ci:‘ep;ancy  include  the  possibility  that  the  actual  hydrograph  for  tbs  event 
may  not  be  entirely  accurate;  and  the  possibility  that  certain  tables  and 


graphs  utilized  by  the  SCS  method  may  not  reflect  the  proper  values  which 
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should  be  used  in  I-bntanae 

3  “  Successful  u^e  of  the  6C3  ir.etlica  rec;uire3  careful  definition  cf 
v.’atershed  parameters  and  flood  routing  ccaff icients .  V/hen  a  ratisfactor 
model  is  finally  obtained,  the  method  represents  a  very  effective,  rapid 


method  for  deterr.inirg  T^eak  flow  rates  to  be  exT>ected.  from  an  actual  or 


■  hji)Dthetical  storm 
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Chapter  X 


KXCOKr«-£XDATIOrIG  AND  SUI-AARY 


bince  xne  nicaeJ.  oi 


Creek  v^aliernbeci  has  now  been  const  re  acted 


and  partially  adjusted,,  further  adjustment  should  be  made  an  more  run-- 
off  even  As  occur  on  Ixick  Creek,  Also^  when  adclitional  discharge  measure¬ 
ments  are  made  by  the  USGS;  the  reliability  of  the  present  rating  curve 
at  the  Idontana  State  University  recorder  on  Duck  Creek  should  be  checked, 
M)del£  of  other  Montana  watersheds  should,  be  constructed  to  ascer¬ 
tain  the  variability  in  runoff  characteristics  anoxig  theiric  Duck  Creek 
may  not  be  representative  of  eastern  Montana  w’atersheds  for  soiue  uni’ ore- 
seen  reasons,  Tlius,^  v’ith  other  iacdG3.s  avrolable^  the  errors  in  design 
storra  peak  disclvarge  predict  ions  would  be  diminished. 

Another  area  of  potential  study  W’ould  be  the  tooting  of  the  various 
facets  of  the  SCS  inethod.  Since  the  present  progran  is  based  largely  on 
the  method  used  before  the  advent  of  the  digital  cornputor  program;  there 
are  components  of  the  method  that  could,  be  more  sophisticated  without 
making  the  problem  solution  more  difficult.  It  would  appear  that  the  syn 
thetic  hydrograTJh  s’nape  could  be  irnde  more  flexible  to  account  for  dif¬ 
ferences  in  area;  shape  of  X’atershcd;  climatic  condition  and  locatj.on, 
Calcu.lation  of  the  theoretical  time  of  concentration  and  lag  tine  needs 
further  frtudy.  The  iv.rportance  of  synthesiaing  the  timing  of  runoff  has 
been  shown  by  the  exa Tuple  used  herein. 

In  order  to  make  further  study  both  easier  and  more  profitable;  the 
SCS“PF}?  should  be  altered  s3.ignt3.y  for  use  on  the  new  digital  computer 


presently  bej.ng  Instailed  at  Montana  State  University, 
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SuKnary 


Prediction  of  the  peak  discharge  of  a  rain-caused  runoff  event  frosi 


a  watershed  by  the  'jiethod  developed -  by  the  SC3  'has  been  discussed. 


u.  l.U 


development  of  this  rictbcl  uas  outlined  in  Chapter  IIIo 


liore  recently  the  ructhoa  has  beez'i  progrc?.riovcd  for  use  vith  an  elec 


tronic  computer.,  A  description  of  this  program  and  an  ezq.hlanaGion  of  itts 
operation  is  found  in  Chapter  IV. 

Ibe  ability  of  the  SCS  method  to  accurately  predict  'peak  discharge 
was  -partially  tested  with  use  of  the  computer,  line  model  of  the  Duck 


Creek  watershed  was  constructed  and  v’as  subjected  to  the  storm  of  June  l6 


1965.  Although  tlia  test  run  is  noi;  sufficieru}  to  evaluate  the  prelecting 


ability  of  the  program^,  the  model  Is  now  available  for  Duck  Creek  e.nd  can 


be  further  adjusted  by  testing  other  storms 


APPENDIX 
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APPENDIX  A 

Solution  of  Runoff  Equation 


sor-'OiO'^fOcu  —  o 


Direct  Runoff  (Q)  in  Inches 


Oi 


Rainfall  (P)  in  Inche 


«  Mt  f't  ■  c»  I#  ^ 
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APPENDIX  B 


Hydrologic  Soil  Types 

A.  (Low  runoff  potential).  Soils  having  high  infiltration  rates 
even  when  thoroughly  wetted  and  consisting  chiefly  of  deep,  well 
to  excessively  drained  sands  or  gravels.  These  soils  have  a  high 
rate  of  water  transmission. 

B.  Soils  having  moderate  infiltration  rates  whan  thoroughly  wetted  and 
consisting  chiefly  of  moderately  deep  to  deep,  moderately  well  to 
well  drained  soils  with  moderately  fine  to  moderately  coarse  tex¬ 
tures,  These  soils  have  a  moderate  rate  of  water  transmission. 

C.  Soils  having  slow  infiltration  rates  when  thoroughly  wetted  and  con 
sisting  chiefly  of  soils  with  a  layer  that  impedes  downward  move¬ 
ment  of  water,  or  soils  with  moderately  fine  to  fine  texture.  Thes 
soils  have  a  slow  rate  of  water  transmission. 

D.  (High  runoff  potential).  Soils  having  very  slow  infiltration  rates 
when  thoroughly  wetted  and  consisting  chiefly  of  clay  soils  with  a 
high  swelling  potential,  soils  with  a  pernanent  high  water  table, 
soils  with  a  claypan  or  clay  layer  at  or  near  the  surface,  and 
shallow  soils  over  nearly  impervious  material.  These  soils  have  a 
very  slow  rate  of  water  transmission. 
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APPEI^DIX  C 

Runoff  curve  numbers  for  hydrologic  soil-cover  complexes 

(Antecedent  moisture  condition  II,  and  =  0,2  S) 


Cover 

land  Use  Treatment  Hydrologic  Hydrologic  soil  group 


or  practice  condition 

A 

B 

c 

D 

Fallow 

Straight  row 

— 

77 

86 

91 

94 

Row  crops 

U 

Poor 

72 

81 

88 

91 

Good 

67 

78 

85 

89 

Contoured 

Poor 

70 

79 

84 

88 

II 

Good 

65 

75 

82 

86 

Contoured  &  Terraced 

Poor 

66 

74 

80 

82 

It 

It 

Good 

62 

71 

78 

81 

Small 

Straight  row 

Poor 

65 

76 

84 

88 

grain 

Good 

63 

75 

83 

87 

Contoured 

Poor 

63 

74 

82 

85 

Good 

61 

73 

81 

84 

Contoured  &  Terraced 

Poor 

61 

72 

79 

82 

Good 

59 

70 

78 

81 

Close-seeded 

Straight  row 

Poor 

66 

77 

85 

89 

legumes  l/ 

It  It 

Good 

56 

72 

81 

85 

or 

Contoured 

Poor 

64 

75 

83 

85 

rotation 

It 

Good 

55 

69 

78 

83 

l^feadow  Contoured  &  Terraced 

Poor 

63 

73 

80 

83 

n 

If 

Good 

51 

67 

76 

80 

Pasture 

Poor 

68 

79 

86 

89 

or  range 

Fair 

49 

69 

79 

84 

Good 

39 

61 

74 

80 

Contoured 

Poor 

47 

67 

81 

88 

11 

Fair 

25 

59 

75 

83 

It 

Good 

6 

35 

70 

79 

Meadow 

Good 

30 

58 

71 

78 

Woods 

Poor 

^5 

66 

77 

83 

Fair 

36 

60 

73 

79 

Good 

25 

55 

70 

77 

Farmsteads 

59 

74 

82 

86 

Roads  (dirii ) 

2/ 

72 

82 

87 

89 

(hard  surface  )  ^ 

---- 

74 

84 

90 

92 

1/  Close-drilled  or  broadcast. 
2/  Including  right-of-way. 
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APPENDIX  D 


MANNING  Q 

CALCULATION  FROM  X-SECT I 

r^^  1 

V-'  1  1 

NOTES 

s. 

C 

A 

AREA  BETWEEN 

TWO  X-SECTION  POINTS 

Vw 

AX3EC 

= 

TOTAL  AREA  OF 

ENTIRE  X-SECTION 

C 

ATOTL 

TOTAL  AREA  OF 

PORTION  OF  X-SECTION 

FOR  EXAMPLE, 

AREA  OF  LEFT  FLOOD,  CH 

A  \  H  l' >1 

EL 

, 

OR  RIGHT  FLOOD  PLAIN 

c 

COEFN 

= 

MANNING  N  C 

OEFFICIENT 

/* 

DELWS 

= 

DISTANCE  Y  TO 

WATER  SURFACE  FROM  NEXT 

PO 

IN 

T 

1 

r- 

(X,Y)  ABOVE  V,' 

ATER  SURFACE 

IHWSEL 

rilGHEST  WATER 

surface  to  be  invest  IGA 

TED 

r- 

LCWSEL 

= 

LOV/EST  'WATER 

SURFACE  TO  BE  INVEST  1  GAT 

ED 

NOW 5 EL 

— 

MO.  0-  WATER 

SURFACE  ELEVATIONS  TESTE 

r' 

nOXGA 

- 

NO.  OF  X-SECT 

ION  AREAS 

C 

iXOXSP 

NO.  OF  X-SECTION  POINTS,  10  IS  THE  MA 

XI  MUM 

r 

s. 

XSU3XS 

NEW  SUB  X~SEC 

TION 

c 

IF  - 

0  NEW  SUB-SECT  lOf-i  TO  B 

C  O  F 

1  \  L_ 

D 

c 

IF  = 

1  NEW  X-SECTION  TO  BE 

n  O’  A 

D 

c 

P 

PERIMETER 

- 

r' 

PTOTL 

= 

TOTAL  PERIMET 

ER  IN  SUB  X-SECTION 

C 

PXSEC 

- 

PERIMETER  OF 

WHOLE  X-SECTION 

r 

SAVXL 

= 

LOCATION  TO  5 

AVE  VALUE  OF  X  ON  LE^T  SIDE 

OF  WS 

C 

INTERSECTION 

TO  BE  'IS ED  WITH  flE'-/  'SS  E 

LEV 

AT 

IONS 

c 

SAVXR 

= 

SA.'IE  AS  SAVvL 

EXCEPT  FOR  RIGHT  SIDE 

c 

SLOPE 

CHANNEL  SLOPE 

c 

n 

= 

TOTAL  DISChlARGE  IN  SUB  X-SECT  ION 

c 

QXSEC 

DISCHARGE  FRO 

M  W!!OLE  X-SECTION 

c 

VXSEC 

AVERAGE  VELOC 

ITY  FOR  WHOLE  X-SECT  I  Ofl 

c 

WSELEV 

WATER  SURFACE 

ELEVATION 

c 

X(  I  ) 

DISTANCE  FROM 

HORZ  CONTROL 

c 

Y  (  I  ) 

ELEV.  ABOVE  DATUM 
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c 

c 

c 

c 

c 

c 

c 

c 

c 
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r' 

V* 
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tv- 
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D  C 

L-i  w. 


DISCHARGE  Q  CAN 
IRREGULAR  NATURAL 
MAY  3E  SUBDIVIDED 


calculated 


FR 


OM  MANNINGS 


T  1 

1  rv  i 


BREAK: 
SEC 


3 


4 


THERE  ARE 


O 

INTO 
FOUR 


CR 


0  55- S 
SUB-SECT  I 
DAT/'' 


T IONS. 


r  Z'  c, 

r\  ^  r\  \  ^  O 


h\\\  \  ' 

th: 

NS  AT  SHARP 
NECESSARY 


:quation  for 

■  V/HOLE  X-SEC 

geometrical 

for  EACH  SUB- 


T  n  ^ ! 

i  j  I  ■i  * 


CARD 

1 


I  I 


i-J-RE  FORMAT  IS  AS  FOLLP'-'^ 


1 

6 

■>  1 
1  X 

16 

21 

31 

1 

X 


NS 

80 

5 

10 
15 
20 
3  0 
40 

n 

O 


9-16 
17  -  24 
ETC. 
1  -  8 
ETC. 


^  « 

I ''4  5 


PRINT  ANY 
LOCATIO 
NOXSP 
LOV/SEL 
IHWSEL 
NSUBXS 
COEFN 
SLOPE 


n 


X 


X 


DIET 


IN 

IN 

TO 


TITLE  SHO'NING 
SECTION  NUMBER 
JUST  I  PI  ED 
JUSTIFIED 
JUSTIFIED 
JUSTIFIED 

X  X  X  X  .  X  X  X  X  X 
XXXX.  XX^'XX 


'ATERSH5 


lA  - 

TC. 


D  NAM 


RIGHT 
RIGHT 
IT 


R  I  ( 


IT 

A 


T  ■ 

\  i  VJi 


FORMAT 
FORMAT 
1ST  POI 


^  t  • 


FROM 


REFERENCE 


WITH  FORMAT  XXXXX.XX 


DIST 


T 


END 

.3RD 


PT 


r  w  ■  I 


REFERENCE 


PT 


ELEVATION  Y 


CORRESPOND  I 


IG  TO.  X 


OF 


CARD  3 


AFTER 
ANOTHER 
CPERATI 
FOR  THE 


ALL  5UB-S 


J  J  .O  C 


X' 


c  r  T  T  o 


r 


ARD 


^  1 1  »  ^  T 

‘*lu  o  ! 


r\  \  I 
^  l\  • 


X- 


i  s  H  A  V  E  R  L^N  FOR  G  I  V  c  N 
READ  TO  INITIALIXE  THE  SUI-IMARY 
THE  CARD  MUST  HAVE  THE  LOWEST  VALUE  OF 
SECTION  IN  COLUMN  1-5,  RIGHT  JUSTIFIED. 


n  r- 

D 


PRINT 

.OWSEL 


V 


I 


i 


A 


t| 


I 

I 


i  u 


-  - 


.  .  .« 


'  ■  t*  „ 

’  T'  ■ 


t'^ 

r,  J 

1  /  ■  t 

:>2^ 

-X  ' 

j* 

:  r  .-: 
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_  I 
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ip* 


^  r 

A  i 
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..M 

‘  .  *  ^  . 

TJff  -r  r^i-g  X 


-  ii 


‘1  ..i’'t.-^ 

I 


C  •  «i«  1  I* 


'  t 

m  W 


^  % 
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r.  cd  T 
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nv 


t  ♦  4-# 


):^r 


4 

I* 


-  c 


> 

:) 

D 

:> 


4 

?» 


if 


“96- 


D  I  MF:^5  ion  T  I  TLE  {  4  3  )  ,  X  (  1  3  )  ,  Y  (  1  ^  ,  DXD  Y  (  9  )  ,  A  (  9  )  ,  P  [  9  )  ?  D  (  1  3  ) 
1  VX5EC  (  6''  "  )  ?  AXSrc  (  63'^  )  ,  PYSFC  (  603  )  ^OXOFC  {  6'^0  ) 

PRINT  104 

104  rCR‘'AT  {  60[i0  '  NANNING  Q  CALCULATION 


1X~SECTI0N 


M  7  rr  c 

1  ‘4  1  L_  J> 


c 

c 

c 

c 


NOTE  THAT  ThIE  FOLLOWING  DO  LOOP 


15  GOOD 


ONLY  FOR 


TlOf^S  FROM  2500  TO  3100  FEET  CHANGE  ELEV.  Ill  STA 
19  TO  FIT  'WATERSHED,  DIMENSION 


CHANGED  FOR  VXSE'^ 


AX3EC,  PXSEC,  AND  OXSEC. 


•c.  n  /  s  I  '4  o  L_  L_  L_  V'  •  i  I  j  o 

STATEMENT  MUST  ALSO 
M  j  /.  .  C  O  N 


FROM 

ELEVA- 

TEMENT 

3E 

STARTS 


19 
2  4 


89 


IN  STATEMENTS  24,  106,  ANC 
DO  23  KK=2501,3130 
K-:<K~2  5  30 
AXSEC { N 


107  MUST  CONFORM 


N  >  ^  »  - 

,  _  'A 


9! 


Yi 


PX5EC{K)=^3. 

QX5EC ( K ) =0 .3 

READ  89, ( T  ITLE(  I ),  1  =  1 ,43) 

FORMAT ( 40A2 ) 

PRINT  90  ,  ( T  ITLE (  I  )  ,  1  =  1 ,40  ) 

FORMAT ( IHl ,40A2/// ) 

READ  91 ,NOXSP  ,LOWSEL , IHWSEL , NSUBXS ,COE^N ,SL^^2 


FORMAT ( 41 5 ,2F1 
I  F { NO XSP 


n 


2  READ 


92 


16,16,2 

,  ( X {  I  )  ,  1= 1  ,NOXSP) 

C  D  1 


READ  92 , ( Y(  I )  , I- 1  ,NOXSH 


O  9 


93 


10  3 


1  =  1  ,NOXSP ) 


FORMAT { 1 :F3.2 ) 

PRINT  93 ,  ( X{  I  )  , 

FORMAT  (4H  X=  l'^FS.2/) 
PRINT  103 ,  ( Y(  I)  ,  1  =  1  ,N 
FORMAT (4M  Y=  13F0.2//) 
NOX5A=NOXSP“ 1 
PRIiNT 


OX. I 


P  ) 


r>  n 

7  , 


r  o  r  p  ,M 
V  u  i_  r  M 


9  5  FOR  M  2\  i  (  9  H  3 


I'i 


^PE 


PRINT  97,SLO 
97  FORMAT!  9H  SLOPE 
PRINT  101 

101  FQR'V,T(62H  WATER 

1  VELOCITY  PERIMETER 

2  (SO  FT)  (EPS) 

DO  15  J  =LOWSEL,  IH'-;SEL 

ATOTL-3 .0 

PTOTL=0. 3 

DO  3  I=l,NOXSP 

AJ  =  J 

D( I ) =Aj- 
I  F  (  D  (  I  ) 

4  D (  I  } =0 
3  CONTINUE 


=F10. 5/ ) 

=  F 1  3 . 5  /  /  /  ) 


SURFACE  DIE 


CHARGE  Q  EN8 
/61H  •  ELEVATION 

( FT )  // ) 


AR 


■  A 
.  M 


i  K  X  J 


4 , 3 


(  CFS) 


I 


' 

: 


( 


-  97  - 


r- 


C 


i  / 


DO  18  I  =  1,00X5 A 
I  I  -  I  +  1 

IF(  Y(  I  I  )~-Y(  I  )  )  5,17,5 
DXDY(I)=C 
GO  TO  16 

5  DXDY (  I  )  =  ( X (  I  I  )-X{  I )  ) / ( Y(  I  I  ) -Y (  I  )  ) 

18  COMTIO'JE 

CALC'JLATE  X  AND  Y  COORDINATES  WHERE 
X-5ECTI0N  0;^  LEFT  SIDE  '"AGING  DOWNSTREAN 


INTERSECTS 
THE  IF{ 


BYPASSES  THIS  SIDE  IF  X-SECTION  IS  0^  RIGHT  FLOOD  S 
TlOr^  ONLY 


t/  ^  1 
t  \  X 


IF(DXDY(  :<)  )  6,31,31 
81  SAVXL  =  X(;<) 


KK=;< 

GO  TO  9 

6  K=L-+  1 

IF  (  D  (  :<)  )  6,6,7 

7  KA-:<-l 
SAVXL  =  X(K'<) 

DELWS=Y{KK)~AJ 
IF  (DELWS)  9,9,8 

8  X{KK)=X(KN)-DELWS-DXDY(KK) 

CALCULATE  X  AND  Y  COORDINATES  WHERE  .  S .  I'lTERECT 
X-SECTICN  ON  RIGHT  SIDE  FACING  DCWNSTREAN 

9  L=NCXSP 
LLL=L-1 

I F ( DXDY( LLL ) )  82,82,10 
32  SAANXR^XC-.) 

LL  -L 
L  =  LLL 
GO 

10  L  =  L-~1 


TO  13 


r;  r\ 


IF  ( D ( L  )  )  10,10,11 

11  LL=L+1 
SAVXR-X( LL ) 

DELV'5  =  Y(  LL)-AJ 
IF  (DELWS)  13,13,12 

12  X( LL ) =X ( LL ) -DELWS-DXDY ( L ) 
calculate  AREA  ,  pern  IMETER  ,  AND  ^AANNING 

13  DO  14  I=KK,L 
11=1+1 

A  (  I  )  =  (  AD  S  F  (  X  (  I  I)  -  X  (  I  )  )  }  -'M  (  D  (  I  I  )  +  2'  (  I  1  !  /  2  .  ) 

P(  I  )  =SGRTF  (  (  (  A3SF(  X(  I  I  )  -  X(  I  )  )  )  ---^'-2  .  )  =  (  (  ADS"'  (  D(  I  )  -D  (  I  I  ) 
l-«--"2  .  )  ) 

ATOTL  =  ATOTL  +  A (  I  ) 

14  PTOTL  =  PTOTL+P (  I  ) 


o  u.i 


-  98  - 


106 


15 

2  0 
10  5 


C 


{  1 .45/CCZF,0  {  (  ATOTL-X--X-1  .  667  )  /  (  PTOTL^'--^"0 . 667  )  )  "  (  SLOPE 

1  "Ir  0*5) 

AVEVEL=Q/ATOTL 
X(LL)=SAVXR 
X{ KK) =SAVXL 

PRIXT  100,  J,0,AT0TL,7VEVlL  ,PTCTL 

FORMAT  (  3X  ,  15  r  8X  5  no.  ^  ,  3X  ,  1  ?  3X  ,  ^6. 2  »  6X  ,  ^6 . 1  ) 

L-J-2500 

PXSEC ( L ) -PXSEC { L ) +PTOTL 
AX5EC ( L ) =AXSEC ( L ) +ATCTL 
QXSEC ( L ) =QX5EC ( L ) +Q 
VXSEC ( L ) =0X5EC ( L ) /AX5EC ( L ) 

IF(N5UBXS)2C ,1,20 
PR  I, XT  105 

F0RMAT(3oRl  SUMLiARY  FOR  EXTIRE  CROSS  SECTIO?>l  ///) 


PRI  X 

T 

1 

1 

READ 

22 

1 

U 

r\ 

ViS 

r 

I 

2  2 

FORM 

A 

T  ( 

I 

6 

) 

DO  2 

1 

J 

= 

I 

L- 

0  l'/ 

L_ 

L  , 

IH 

1 ,  f  r 

"l 
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L-J- 

2 

0  J 

-J 

21 

PR  I  X 

T 

1 

1 

XS 

LI 's- 

{  L  ) 

. AXSEC { L )  ,  VXSEC 

GO  TO 

1 

> 

PL  AC 

r" 

r_ 

T 

\ .  * 

o 

lJ 

LA 

fC 

n 

-7  A'\ 

TA 

^  A  r>r^  c  at  [rrir\  otr 

'  V  /■■'  1  _  1  <  Lx  '  J  . 

16 

CALL 

EX 

r 

T 

O 

»  r  ,  \ 


S  E  C  (  L  ) 


l-N  A  T 

w’  *-1  1  /  \ 
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